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Abstract 

We compute analytically the tree-level annihilation rates of a collection of non- 
relativistic neutralino and chargino two-particle states in the general MSSM, in- 
cluding the previously unknown off-diagonal rates. The results are prerequisites 
to the calculation of the Sommerfeld enhancement in the MSSM, which will be 
presented in subsequent work. They can also be used to obtain concise analytic 
expressions for MSSM dark matter pair annihilation in the present Universe for a 
large number of exclusive two-particle final states. 



1 Introduction 



The presence of dark matter (DM) in the Universe is one of the few empirical evidences 
that the current Standard Model of particle physics cannot be complete. The dark 
matter density ^dm = 0.111(6)//i^ {h = 0.710(25)) P is now determined very precisely 
from various observations. While the particle nature and genesis of the dark matter relic 
density remains unknown, it is intriguing that it can be explained naturally through 
thermal production and freeze-out of a particle with electroweak interaction strengths 
and mass of order of the TeV scale. In this simple scenario freeze-out occurs when 
the Universe cools below the mass of the particle and the DM particles become non- 
relativistic, with typical velocities of order v ~ 0.2 c. The DM pair-annihilation cross 
section, which determines the relic density, can then be expanded in a Taylor series in 
V, and keeping only the first two terms is usually a very good approximation: 

O-ann^^rcl ^ a + b V^^^ . (l) 

Here v^e\ = \vi — V2\ denotes the relative velocity of the annihilating particles in their 
center-of-mass frame. Furthermore, when dark matter particles pair-annihilate in the 
present Universe, potentially revealing themselves in cosmic ray signatures, the typical 
velocities are v ~ 10~^c, and the annihilation occurs even deeper in the non-relativistic 
regime. 

Among the many models that contain weakly interacting dark matter candidates at 
the TeV scale, the minimal supersymmetric standard model (MSSM) has been studied 
most extensively, and for quite some time [213] . Several programs are available [HJl] which 
compute the annihilation cross section of the lightest neutralino, together with possible 
co-annihilation processes, in the MSSM numerically in the tree-level approximation. The 
observed relic density then provides a valuable constraint on the parameter space of the 
model, complementary to those from collider physics. Given the precision of Qbm, it 
seems desirable to compute the cross section parameters a, b at the one-loop level. This, 
however, is a daunting task due to the complexity of complete one-loop calculations in the 
MSSM, and the number of individual annihilation processes that add up to the total cross 
section. Nevertheless, such calculations have been performed for certain scenarios where 
QCD corrections are the most important ones [6]-E], or in certain approximations [9|[T0]. 
The calculation of the full electroweak corrections has been started [TTlfT^ . 

There exist situations when quantum corrections become exceedingly large and can- 
not be neglected. In non-relativistic scattering and annihilation of DM particles this 
happens when the Coulomb (Yukawa) force generated by massless (massive) particle 
exchange between the DM particles becomes strong at small relative velocities, a phe- 
nomenon also known as "Sommerfeld effect". In the MSSM this situation is naturally 
realized when the lightest neutralino (LSP) has mass above one TeV, in which case the 
neutralino is almost a pure gauge eigenstate, and degeneracies and co- annihilation effects 
in the neutralino-chargino sector are generic. The Sommerfeld effect in the MSSM was 
first studied in the wino- and Higgsino-limit by Hisano et al. [T3l[T¥j , and subsequently in 
"minimal dark matter models" [15] that resemble the MSSM in the above-mentioned lim- 
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its. In these heavy dark matter scenarios the annihilation cross section can be enhanced 
by more than an order of magnitude, since the typical distance of DM particles at small 
velocity is within the long-range part of the Yukawa potential generated by exchange 
of the electroweak W and Z gauge bosons. The suggestion [16] that the Sommerfeld 
enhancement due to the exchange of a new, light particle may generate an excess in the 
cosmic ray positron spectrum has generated a resurge of interest in this effect. 

The present work aims at improving the calculation of the dark matter annihilation 
cross section and relic abundance by including the Sommerfeld radiative corrections in 
the general MSSM, beyond the previously considered wino- and Higgsino-limit. The idea 
is that even when the Sommerfeld correction is not of order one, unlike in scenarios of 
TeV scale LSPs, it may still constitute the dominant radiative correction in a significantly 
larger portion of the MSSM parameter space. Related work has been undertaken recently 
in [T7H20] . Our approach differs from or extends these works in several aspects. 

• We use the non-relativistic effective theory approach to separate the short- distance 
annihilation process from the long-distance Sommerfeld effect, which is encoded in 
the matrix elements of local four-fermion operators. The approach is very similar 
to the NRQCD treatment of quarkonium annihilation [21], except that we deal with 
scattering states of several species of particles interacting through the electroweak 
Yukawa force. 

• Since electroweak gauge boson exchange may change the two-particle state (for 
instance, scatter a neutralino pair into a pair of oppositely charged charginos), the 
annihilation process is described by a matrix in the space of two-particle states, 
which is not diagonal. The off-diagonal terms cannot be obtained from the tree- 
level cross sections computed by numerical programs, and have not been considered 
previously, except in the simplified situation of the strict wino- and Higgsino- 
limit [I3HIS1I2U]. 

• We compute the expansion of the short- distance annihilation cross section ana- 
lytically rather than numerically. The only systematic previous analytic calcula- 
tion [22] refers to the annihilation of two LSPs, but does not include co-annihilation 
channels and the above-mentioned off-diagonal annihilation matrix entries. 

• The non-relativistic annihilation cross section can be organized in a partial wave 
expansion. The leading-order term a in ([T]) contains the leading-order contributions 
from S-wave annihilations, whereas the second term b encodes both S- and P-wave 
annihilation contributions. The Sommerfeld correction factor is different for the 
S- and P-wave contribution. Its consistent implementation therefore requires a 
separation of b into its two components, which has not been done before, but can 
be implemented relatively easily within our analytic framework. 

The present paper is devoted to the analytic calculation of the dominant S'-wave annihi- 
lation coefficient a and to outlining the general framework. The subleading term b and 
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the calculation of the Sommerfeld effect in the MSSM with almost degenerate neutralinos 
and charginos will be presented in [23] and [J^, respectively. 

The outline of this paper is as follows. In Sec. [2] we introduce the effective La- 
grangian framework and the method of calculation. In the non-relativistic MSSM the 
short-distance annihilation process is encoded in the Wilson coefficient of a four-fermion 
operator. We introduce the required notation and discuss the expansion in the mass 
differences of the nearly degenerate neutralino and chargino species. Sec. [3] discusses 
various examples of annihilation cross sections obtained from our analytic calculation. 
In this section we also perform checks by comparing some diagonal annihilation matrix 
entries with numerical cross sections and with [22]. The complete analytic results are 
rather lengthy. We provide them in appendix El together with the rules to construct the 
MSSM coupling factors of the various diagrams. In Sec. IHwe explain why it is convenient 
to employ Feynman gauge despite the fact that this requires the computation of many 
unphysical final states. We also illustrate the importance of including the off-diagonal 
annihilation matrix entries in the computation of the Sommerfeld-corrected cross section 
on the example of a heavy wino-like MSSM parameter point. We summarize in Sec. O 



2 Effective Lagrangian and method of calculation 

2.1 The Lagrangian in the effective theory 

We describe the kinematics and interactions of neutralinos and charginos moving at small 
velocities within a non-relativistic effective theory (EFT), the non-relativistic MSSM 
(NRMSSM), that contains only nearly on-shell non-relativistic chargino and neutralino 
modes, while the effects from higher mass and virtual modes are encoded in the Wilson 
coefficients of higher-dimensional operators. The neutralinos and charginos described in 
the EFT approach are those whose masses are nearly degenerate with the mass mLSP of 
the lightest neutralino. The corresponding effective Lagrangian is given by 

^NRMSSM ^ ^^.^ ^ ^^^^ ^ ^^^^^ ^ j^igj^g^ Q^^g^ ^g^^g _ (2) 

The kinetic part of the Lagrangian for no < 4 non-relativistic neutralino species and 
< 2 non-relativistic chargino species is given by 

no 
i=l 




+ 5Z 5Z ^1 I - ^"^^ - "'LSP) + Tj;;;^ ) ipj ■ (3) 



LSP 



The fields C,i and ipj = f]j,Cj represent the non-relativistic two-component spinor fields 
of non-relativistic neutralinos (Xi) and charginos {xj and x/)? respectively. This EFT 
setup with one reference mass scale, mLSP, is suited for the description of (neutralino) 
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dark matter annihilation processes in the present Universe as well as for the computation 
of dark matter co-annihilation reactions with further nearly mass-degenerate neutralinos 
and charginos in the context of the relic abundance calculation. The EFT framework 
can easily be extended to the case where the non-relativistic particle species are (nearly) 
mass-degenerate with respect to two distinct scales m^^f\ with m^^J- ^ '^rcf- that 
case, the mass differences {nik — "n^LSp) in (El) have to be replaced by nik — m^ci,k, where 
each mref,fc is given by one of the scales m).^^ ■ In that way, an entirety of hydrogen-like 
two-particle states can be described, within which a set of light, nearly mass-degenerate 
and another set of heavy, nearly mass-degenerate particles exists. Our results for the 
absorptive part of the Wilson coefficients, specified in Sec. l2.2l and given in the Appendix, 
cover both, the cases of a set of particles nearly mass-degenerate with the neutralino LSP 
and a set of non-relativistic hydrogen-like neutralino and chargino systems. 

The term £pot summarizes (instantaneous) Yukawa- and Coulomb potential inter- 
actions that arise through the exchange of SM gauge bosons and Higgs particles. The 
generic form of £pot reads 

^pot = - j d^f $^;(x, r ) Vijki{r) $ij(x, f ) (4) 

where the fields describe a two-body state of the form Xi'Xj'; 'X^xf xfxf and a sum 
over repeated indices is implicit. Vijki{r) thus represents the potential interactions among 
two-body states (ij) and (kl) in configuration space, with r the spatial 3- vector denoting 
the relative distance in the two-body system. The explicit form of the potentials between 
neutralino and chargino species for a given MSSM point will be given elsewhere |24j . 



2.2 Basis of the dimension 6 operators in 6Cann 

Within the NRMSSM, we aim to describe neutralino and chargino pair annihilation pro- 
cesses into two-particle final states of Standard Model (SM) and (light) Higgs particles, 
which are not non-relativistic. The theory will contain effects from virtual and higher- 
mass Higgs and SUSY particle modes as well, encoded in the EFT operator-coefficients 
and parameters. The specific case of resonant s-channel pair-annihilation reactions can 
be covered by adding a resonance width in the analytic results that we give in the 
Appendix. Yet we exclude the case of accidental mass degeneracies of further SUSY 
particles with the set of non-relativistic neutralinos and charginos. 

The SM and light Higgs particle final states in the neutralino and chargino pair- 
annihilation reactions are not described within the non-relativistic effective theory, as 
they are characterized by velocities outside the non-relativistic regime. However, since 
the hard inclusive pair- annihilation processes take place within distances of order 1/mLSP, 
we can incorporate the short- distance annihilation rates of non-relativistic neutralinos 
and charginos in the effective theory through the absorptive part of Wilson coefficients of 
local four-fermion operators in (5£ann, following the approach of [21]. The full annihilation 
rates in the non-relativistic effective theory are given by the absorptive part of the matrix 
elements of these four-fermion operators. While the matrix elements of the operators 
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Figure 1: Diagrammatic picture for the relation among the annihilation amplitude and 
the absorptive part of the corresponding forward scattering amplitude in presence of 
long-range potential interactions. 



themselves may encode long-distance effects, giving rise to Sommerfeld enhancements, 
the contribution to the hard annihilation reaction factors out in the form of the Wilson 
coefficient. 

In contrast to the application of this formalism to quarkonium annihilation in QCD 
[2Tj . we are going to describe annihilations of scattering states instead of bound states 
and allow for more than one non-relativistic particle species. The latter allows for the 
possibility, that (long-range) potential interactions (indicated by the grey oval in Fig. [1]) 
lead to transitions from the initially incoming two-particle state XiXj to another two 
nearly on-shell non-relativistic two-particle state XeiXe2 prior to the annihilation reaction. 
Unitarity relates the phase space integrated product of annihilation amplitudes XiXj ~^ 
XaXb in the first line of Fig. [T]to the absorptive part of the forward scattering amplitude 
XiXj ~^ XiXj depicted in the second line, where XaXb generically denotes a pair of SM 
and light Higgs particles. Note that due to the presence of the long-range potential 
interactions, the hard annihilation reaction is determined by the absorptive part of the 
XeiXe2 ~^ XeiXez amplitude, as can be seen in the second line of Fig. [TJ It is worth to 
stress, that the XeiXe2 particle pair is not necessarily equal to the Xe^Xe-i pair, such that 
apart from true forward scattering reactions XeiXe2 ~^ XeiXe2! we encounter off-diagonal 
XeiXea ^ Xe^Xes rcactlous as wcU. 

In this paper we are concerned with the calculation of the absorptive part of XeiXe2 ~^ 
XaXb — )■ XciXe-j, 1-loop reactions, encoding the hard tree- level XeiXe2 and Xe4,Xe3 anni- 
hilation processes. The full annihilation rates, including the long-range effects shall be 
studied elsewhere [23]. To determine the absorptive part of the Wilson coefficients, we 
evaluate the absorptive part of the hard XeiXe2 ~^ Xe^Xe-j, 1-loop scattering amplitude 
within the MSSM and match the result with the tree-level matrix element of four-fermion 
operators contained in 5£ann in the effective theory. At tree-level the annihilation rates 
can be given separately for every final state XaXb-, since the tree- level processes are 
free from infrared divergences. In higher-orders the formalism applies to the inclusive 
annihilation cross section [21], or to suitably defined infrared-safe final states. 
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neutral reactions 


single-charged reactions 


double-charged reactions 




x^x^ X^X^ 
x'x^ 


x^x^ 

x~x~ x~x~ 



Table 1: Collection of all XexXe2 Xe4,Xez scattering reactions. The labels Cj on the 
fields Xei are suppressed in the above table. If Xe^ represents a field Xg . , the label e, can 
range over = 1, . . . , ng, whereas = 1, . . . , n+ for the case of a field. 

The leading-order contributions in 5£ann are given by dimension-6 four-fermion opera- 
tors. For instance, the specific dimension-6 four-fermion operator that encodes scattering 
of a non-relativistic incoming neutralino pair XiXi iii a ^5*0 partial- wave state into an 
outgoing XiXi state in the same ^5*0 partial- wave configuration, is given by 

''''-ann ^ ^ •^{11}{11} ^ '^0'' ^1^1 ^1 ^1' 

where the spinor is the charge conjugate of ^, = — ia^^*, and specifies the 
second Pauli matrix. Note that represents the Lorentz invariant bilinear built from 
the non-relativistic particle field which destroys the incoming state of two identical 
X? particles. The factor 1/4 denotes a normalization factor which compensates the 
symmetry factors arising from the number of identical contractions in the tree-level 
XiXi XiXi matrix element. The symbol ffii^^^ ^ (^Sq) denotes the Wilson coefficient 
corresponding to the dimension-6 operator. We can generalize the above expression to 
include all possible spin-0 and spin-1 S'-wave four-fermion operators at leading order 
in the non-relativistic expansion in. Written in a compact form, the contribution of 
dimension-6 operators in 5£ann reads 

CA.d=6 _ \ fXX-^XX (28+1 q \ c^XX^XX (2s+lq \ (p\ 

ami / J / J ^ •'{6162 }{e4e3} ^ '^■^} '^{e4e3}{e2ei} I ) ' 

xx-^xx s=0,l 

where J = s for the case of S'-wave operators considered here. The first sum, taken over 
all non-relativistic 2 — 2 neutralino and chargino scattering processes xx ~^ XXi implies 
the consideration of neutral scattering reactions as well as single- charged and double- 
charged processes. The xx ~^ XX reactions that we take into account are summarized in 
Tab. dJ The spin of the incoming and outgoing two-particle states can be either s = or 
s = 1, such that the terms in the above Lagrangian describe ^5*0 and ^5*1 partial- 

wave scattering reactions. The ff^^^^ i'^^'^^Sj) denote the Wilson coefficients, that 

a •' 16162116463} \ ■J y ' 

correspond to the four-fermion operators ^{e4e^{e2ei} (^^^^'^'j). The indices ei and 62 (63 
and 64) refer to the neutralino or chargino species of the incoming (outgoing) particles, 
and take the values 1 to uq for neutralino species and 1 to for chargino species. 
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XeiXe2 ^ Xe4Xe3 



x° x° ^ x° x° 

x° x° x'x'^ 
x~x^ x° x° 



x^x"* 



x°x^ 



X X° ^ X x° 



X X X X 

x^x" ^ x^x'" 



/nxx^xx (ic 



Se4 Sea 

„t re 

'lei Ses 
Se4 Sea 
'£4 Sea 



Se2 Sei 
Se2 Sei 
Se2 '/ei 
Se2 '/ei 



ft /-c /-ct f 
Se4 Sea ^62 Sei 



se4 



^ Cl ^ei 



'£4 ^ea 



At AC Act A 
Se4 Sea Se2 Sei 



n^X^XX (Z q ^ 
'^{e4e3}{e2ei} '^l^ 



ft ffCc fct^f 

Se4^Sea Se2 Sei 

t ^^c fct^f 

/e4^ Sea Se2 Sei 

Se4"Sea Se2^ '/ei 

^l^^Cea Cel^^ei 



el4^Ce=3 



Se2" Sei 



fct 
Se2 



Se4 Sea Se2 Sei 
^t^^ea ^e1'?^ei 



Table 2: Four-fermion operators for leading-order S-wave XeiXe2 ~^ Xe^Xez transitions. 
The indices e^, i = 1, . . . , 4 on the x-fields are suppressed in the first column. In addition 
to the specified operators there are redundant ones, which are obtained by interchanging 
the field-operator symbols ^,1] or ( (but not the labels) at the first and second and/or the 
third and fourth position in the operator C)xx^xx example, for ^5*0 x'^X^ ~^ x'^X^ 
scattering one of the three field-interchanged operators reads (l^ ^'^^ Cer 



Note that the order of the labels Cj on the Wilson coefficients and the operators is not 
accidental in (|6]). The labels on the operators are given in the order, in which the 
field operators with label occur in the operator. In case of the corresponding Wilson 
coefficients, the indices refer to the actual scattering reaction XeiXe2 Xe4Xe3 5 that is 
described by the operators. For the basis of the operators see Tab. O The x iii the labels 
XX ~^ XX of the operators and Wilson coefficients in ([6]) should indicate the particular 
particle species x° and whose XeiXe2 ~^ Xe^Xez scattering reaction is described, see 
Tab. [TJ A summation over the indices Cj is implicit in (|6]). The normalization factor 
1/4 in ([6]) ensures that the tree-level transition matrix element for ^S'o-wave scattering 
is given by 

(XiXfcl j ^ - /fe^^}fe4ea}(^'S'o) Cg7a}{e2ei} (^'^'o) W |XiXi)|tree 

XX^XX 

= (2vr)^5(^) fe„ - 2 ff^^^^klCSo) (7) 

for all XiXj ~^ XiXk reactions at leading order in the non-relativistic effective theory. A 
similar relation for the tree-level transition matrix element of ^S'l-wave scattering in the 
effective theory holds for all XiXj ~^ XiXk reactions. Note that in order to derive ([7]) one 
has to take into account relations among Wilson coefficients of different operators, which 
will be deduced in the next paragraph. 
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There are redundancies in SC^^^, ([6]), as several operators can describe one specific 
scattering reaction with a Xei and a Xe2(Xe4 and Xes) particle in the initial (final) state. 
This redundancy is associated with operators that arise from interchanging the single- 
particle field operators at the first and second and/or third and fourth position in a 
given O^^^^^. The corresponding Wilson coefficients are related to each other, as they 
encode the same information on a given specific scattering reaction. Consequently, the 
redundancy manifests itself in symmetry relations among the Wilson coefficients under 
exchange of the labels ei ^ and/or 64 -H- 63. These relations read 

•'{e2ei}{e4e3} I J ) — 'h i {eie2}{e4e3} I "^JJ ' 

/•XeiXe2^>Xe3Xe4 /2s + l Q \ /■Xei Xe2 ^Xe4 Xeg (2s+lq \ /o\ 

•'{eie2}{e3e4} I ^J) - 'Is J {e^ei} {e^ei} V J ) ' \°) 

with 

f 1 fors = 
= { ■ (9) 

1 —1 for s = 1 

To exemplify the origin of the first relation in ([8]), let us consider the terms in 5C^^ that 
account for S'-wave x^X^ ~^ XX reactions at leading order in the non-relativistic velocity 
expansion: 

El /■X^X'^-^XX (2s + lQ \ c^XX fcf f /-|r)\ 

4 -I {eie2}{e4e3} { ) ^{aei} ^£2 , {J-^) 

ei,...,e4 

where the operator C^^^gg} stands for the two-field operator that creates the outgoing 
state, and T^^^ is given by r*^*"'^-' = l2x2 and F*^**^^) = a in case of ^5*0 and ^5*1 operators, 
respectively. Using the definition of the spinor we can write 

Cl r^^) U = Cl cr^r(^) = vs Cl r^^^ ^2 , (n) 

with 77^ defined as in (Q. After renaming the labels ei and 62, the terms in (fTOj) can be 
written as 

Ei„ fxV^XX (2s+lc\ r)XX fct f 
4 VW{e2ei}{e4e3} ^ '^JJ '^{6463} ^£2 • 

ei,...,e4 

Comparing to the original expression in (fTOj) . we arrive at the relation 

fX^X^-^XX (2s + lq \ _^ fX°X°^XX (2s + lq \ (■\cy\ 
•'{e2ei}{e4e3} I J ) ~ {eie2}{e4e3} I J ) ' \^'^) 



This equation as well as the more comprehensive relations in (jS]) imply that the Wilson 
coefficients of ^5*0 operators have to be symmetric under the exchange ei -H- 62, whereas 
Wilson coefficients of ^Si operators are antisymmetric under ei O 62. The same state- 
ment applies to the exchange 63 -H- 64 in case of outgoing states. In processes with 
identical incoming or outgoing particles, the above relations in ([8]) imply the vanishing 
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of the ^5*1 Wilson coefficients. This rephrases the well-known fact that a pair of identical 
spin-1/2 particles cannot build a ^5*1 state. 

Finally, a further property of the Wilson coefficients under the exchange of the particle 
labels is directly inherited from the hermiticity of the non-relativistic Lagrangian: 



•^{eie2}{e4e3} I ^J) 



fXX^XX (2s+l q \ 

/{e4e3}{eie2} I ^J) 



(13) 



2.3 Matching condition 

The Wilson coefficients of the four-fermion operators in 5£ann are determined by the 
matching condition 

-^(XiXj X/Xfc) ImSSM, pcrturbativc = ^ f{^v^2}{e4e3}(^'^~^^^j) 

X {XlXkl ^{e^4^}{e2ei}( ^j) \XiXj) InRMSSM, perturbative • (14) 

For this equation to hold, we have to use the same (non-relativistic) normalization of the 
incoming and outgoing states in both the full theory and the NRMSSM. Here we will 
determine the contributions to the Wilson coefficients that describe the tree-level anni- 
hilation reactions into exclusive SM and light Higgs two-body final states XaXb, which 
we shall denote as f^^^^^^s-^xx^^s+i j^^y rjj^^ unitarity of the S-matrix at the diagram- 
matic level establishes a relation among the tree-level annihilation rate for XiXj ~^ XaXb 
and the imaginary part of the 1-loop forward-scattering reaction XiXj XaXb — > XiXj'- 

j [dPS^fi] \A{x^XJ ^ XaXb)\^ = 2 !3 [A{x^XJ ^ XaXb ^ XiX,)] (15) 
4 

We generalize this and define the absorptive part of amplitude A{xiXj ^a^b XiXk) 
as well as the absorptive part of the Wilson coefficients in the following way: 

j [dPS^B] AixiXj ^ XaXb) X AixiXk ^ XaXbY 

= 2 [A{XiXj XaXb XlXk)] labsorptive 

" ^ 4 AeSniV^^'^^''^^^'^) {XlXk\Ofe^,^^e2ei}C'^^Lj)\XiXj) , (16) 

where we have introduced the notation 

fXX^^A^S-S>XX/'2s+l r \ _ rXX^^A^B^XX/^s+lT \ \ (-,fy\ 
•/{ijUlk} \ ^J> ~ J {ij}{lk} \ ^J> labsorptive • \i- ' ) 

With this definition, the absorptive part of a Wilson coefficient that encodes a XiXj ~^ 
XiXj forward-scattering reaction coincides with its imaginary part. 



fXX^^AXB^XXf^s+l 
{eie2}{e4e3} V ^ 
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We make use of the defining relations to determine the absorptive part of the Wil- 
son coefficients J^^^^a^b^xx from the product of the full-theory tree-level annihilation 
amplitudes integrated over the final state particles' phase-space, as given in the first 
line of f|T6l) . Technically this is achieved by considering all 1-loop scattering amplitudes 
XiXj ~^ ^aXb — XiXk with a specific SM or Higgs particle-pair XaXb in the interme- 
diate state and by applying the Cutkosky rules to the Xa and Xb propagators. The 
resulting expression coincides with the first line of (fT6l) . To determine the absorptive 
part of the Wilson coefficients, the expression has to be expanded in the non-relativistic 
momenta of the external particles as well as in their mass differences and an appropriate 
spin-projection has to be performed. 



2.4 Expansion in mass differences in 5C 



ann 



To simplify the notation, we shall replace the indices (ei, 62, 63, 64) by (1,2,3,4) through- 
out this section. Further, we shall adopt the convention that particles 1 and 4 in the 
reaction X1X2 — ^ X4X3 share the same reference mass scale m, while particles 2 and 3 
have masses closer to the reference scale m. Introducing two distinct mass scales for 
the particle species allows us to consider pair annihilations of two particles with sim- 
ilar mass (m ~ m), but also pair annihilation of a hydrogen-like two-particle system 
where one of the particles is much lighter (though still heavy enough to be considered 
as non-relativistic). According to these assignments, we define 

mi = m — 6m , m2 = m — 6m , 

m^ = m + 6m , m^ = m + 6fn , (18) 



with 



mi-Fm4 _ m2 + m3 
m = , m = , (19) 

such that the mass differences read 

6m = ^^-^\ 6m =^^^^. (20) 

The results for the Wilson coefficients presented in the appendix adopt the definitions 
(ITBH2U]) . If for a given process XiXj ~^ XiXk if turns out that the reverse condition, 
~ i^k ~ ^ and mj ^ mi fn, is more meaningful given the actual values of the 
masses, one can make use of the symmetry properties ([8]) to relate the Wilson coefficients 
for XiXj ~^ XiXk to those of XiXj ~^ XkXu which would then conform to the prescription 
above, i.e. m would be equal to the average of the mass of the particle associated with 
field 1 and the mass of the particle associated with field 4, m = (mj + mk)/2. Note that 
the mass differences 6m and 6m in (l20l) obviously vanish in case of diagonal scattering 
reactions X1X2 X1X2, such that m = mi and m = m2 in that case. 
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The absorptive parts of the Wilson coefficients are obtained by matching amphtudes 
for the process X1X2 XiXs with on-shell external states. This implies that the energy- 
conservation relation in the center-of-mass system, 

V~s = E,{p') + E,{f) = E,{p'') + E,{p'^) , (21) 



with Eiip"^) = a/ + p 2 , is fulfilled. Using ( IT8] - [20|) and M = m + m, the expansion of 
the energy-conservation relation (l2T]) for non-relativistic momenta p"^ and p'"^ reads 

^ = M-6m-6m+^ + ... = M + 6m + 6m+^ + ... , (22) 
^ 2/i 2/i ' ^ ^ 

where /i = mfn/M and terms of order p^/fi^ and {5m/ M x p'^/fi) have been dropped. 
This can be rewritten as 

-»/2 -*2 

^ = |--25m-2(5m + ... . (23) 

From (123|) we see that a consistent expansion which treats both and p'^ as small 
quantities of the same order requires that the mass differences 6m, 6fn are also formally 
considered of order p^//i in the expansion of the amplitudes. Note that an expansion 
in mass differences is only required for the off-diagonal scattering reactions where the 
incoming and outgoing XX states are different, as 6m = 6m = for X1X2 — ^ X1X2 
reactions. 

The amplitude for a generic process X1X2 — ^ ^a^b — ^ X4X3 then depends on the 
hard scales (m, m) and on the small scales {p'^/fi, p'^//i, p- p'/fJ^, 6m, 6m) ~ 0{fiv^), 
where v stands for the relative velocity in the two-particle system. In the following we 
enumerate the steps to obtain the absorptive part of the Wilson coefficients from the 
process X1X2 ^a^b XiXs^ including the subleading 0{v^) terms, which will be 
presented in [23] . 

1. The absorptive part of the 1-loop scattering amplitude X1X2 XaXb — )■ X4X3 
with a SM or Higgs final state XaXb is computed by applying the Cutkosky rules 
to the Xa and Xb propagators. The result is written in terms of the mass scales 
introduced above, and expanded in the small scales retaining terms up to 0{v'^). 

2. To 0{v'^) the result contains scalar products with at most two powers of p and p' . 
For the spin-1 configuration, the scalar products also involve the spin-polarization 
vectors n and n' of the incoming (X1X2) and outgoing (X4X3) states, respectively. 
The generic form of the result for spin-1 incoming and outgoing states reads 

{ Coi^Si) + ci(35i) 6m + C2{^Si) 6m + c^('Si) p^ + c^{^Si)p'^ ] n ■ n' 
+ c^i^Po) {p- n) {p ■ n) + c^{^Pi) [p, nf [p' , nf + Cj('P2) p^' n'^ p'^' n^^ 
+ csCS,, 3Pi) n^- [p', nf + cgCP^, 'S,) [p, rifn'" 
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where we have introduced the notation [a, = e'^^^a^hP and a^"^ = a^V + a^V — 
2 a - b6'''^ /3, corresponding to J = 1 and J = 2 Cartesian tensors, respectively. The 
spin-polarization vector n is introduced by replacing the spinor matrix [^^'^^jjj of an 
incoming two-neutralino state by 4j n-aij. Similar replacements apply to outgoing 
two-particle states and states involving charginos. The coefficients Cj are functions 
of m and fn. The ffist term, cq, gives the leading-order contribution, where all the 
others count as (9(f ^). We have further specified the quantum numbers ^Lj of each 
term, which matches the angular-momentum configuration of the incoming state, 
equal to that of the outgoing state except for the Cg-u terms (the ffist quantum 
number between parentheses refers then to the incoming state, the second to the 
outgoing one). For spin-0 incoming and outgoing states, the result simplifies to 

CoCSo) + CiCSo) Sm + cCSo) Sfn + 03(^^0) + cCSo) + c^{'Pi)v-v' • 

(25) 

We have not considered the possibility of spin-0 to spin-1 transitions between in- 
coming and outgoing states in the hard annihilation process, though the transitions 
^Pi and ^-Po,i ^Sq, ^Pi are also allowed at 0{v) by angular-momentum 
conservation. Such spin-changing transitions in the hard annihilation part of the 
full forward scattering amplitude (see Fig. [1]) will also require spin-changing poten- 
tial interactions in the long-range part, in order to bring the spin of the two-particle 
state after annihilation back to the spin of the (left-most) incoming state. Since 
the non-relativistic spin-changing potentials carry an additional ^-suppression, such 
transitions are only relevant for the calculation of the annihilation rates at 0{g'^v'^). 
At present we ignore 0{v'^) effects that arise from subleading non-Coulomb (non- 
Yukawa) potentials and consider only those from the short- distance annihilation. 
Likewise, the terms cg-u included in ( 12^ imply a change of the orbital angular 
momentum which must be compensated by a potential interaction which is also 
f -suppressed in the non-relativistic limit, and can be ignored for our purposes. 

3. By virtue of the energy-conservation relation (1^21) . we rewrite powers of p ^ and p'"^ 
as 

= ^(p'+p") + ^^(5m + 5m) + ... , 
P" = l{p'+p")-'^{Sm + 6m) + ... , (26) 

such that the coefficients multiplying and p'"^ become equal. This convention 
is adopted in order that the Wilson coefficients of the dimension-8 operators with 
derivatives also have the symmetry property (|T3|) under the exchange of the in- 
coming and outgoing states. 
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4. Finally, the Wilson coefficients are identified by comparing the expanded expression 
for the absorptive part of the amplitude A{xiX2 — > ^a^b XiXs) with the 
amplitude for the same process computed with the dimension 6 and dimension 8 
EFT operators in 5£ann- 

The explicit expressions for the leading-order S'-wave coefficients are provided in the 
appendix. While the Wilson coefficients refer to the inclusive annihilation rates, summed 
over all accessible final states, the calculation is performed for individual final states, 
which are therefore also given separately. Such final-state separated results should be of 
interest to the calculation of primary decay spectra of dark matter annihilation in the 
present Universe. 



3 Results 

We have performed a number of dedicated numeric and analytic checks of our results 
for the absorptive parts of the Wilson coefficients. As these expressions also encode 
the absorptive part of XeiXe2 ~^ ^a^b — ^ XeiXe2 forward scattering reactions, which 
are related to the tree-level annihilation cross section for XeiXe2 XaXb processes 
(see ( !T5|) ). a comparison of the analytic non-relativistic approximation to the tree- level 
annihilation cross section with results from a numeric code can be carried out for the 
diagonal entries of the annihilation coefficients. We discuss this in Sec. 13.11 

In addition, we can relate our analytic expressions for partial-wave separated neu- 
tralino LSP pair-annihilation cross sections to existing analytic results available in the 
literature [T3l[TU|22]. This will be briefly discussed in Sec. 13.21 No checks are available 
in the general case for the off-diagonal entries of the annihilation coefficient matrix. 



3.1 Numerical comparison with MadGraph 

The expansion of the exclusive, spin-averaged center-of-mass frame XeiXe2 ~^ ^aXb 
tree-level pair-annihilation cross section in the non-relativistic momentum p of the Xei 
particles is given by 

^xe,xe2^x^xs ^^^^ ^ f^is^) + 3 f^^s,) (27) 

+ p' {fCPi) + I fCPo) + fCPi) + I fCP2) + gCSo) + I gCSo)) + 0{p') . 

Here Wrei = l^ei — ^^2! is the relative velocity of the XeiXe2 P^-ir and v^^ denotes the velocity 
of particle Xet in the center-of-mass frame of the annihilation reaction. We have sup- 
pressed the superscripts XeiXe2 ~^ ^a^b XeiXe2 oti the Wilson coefficients / in ( 127|1 . 
where these expressions explicitly refer to the exclusive (tree-level) annihilation rates. 
Further note that (1271) contains not only the leading-order S'-wave Wilson coefficients 
f{^^~^^Sj) with spin configuration s = 0, 1 of the incoming two-body system, but also 
includes P-wave and next-to-next-to- leading order S'-wave coefficients (denoted with g). 
For analytic results on those coefficients we refer the reader to [23] . 
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In the non-relativistic limit the relation between the relative velocity v^ei and the 
particle momentum p in the center-of-mass frame of the Xei Xe2 annihilation reaction is 
approximated by 

Vrel = \Ve, ' V,, \ = \p\ f ^ + 0^] . (28) 

Together with fl27|) . this relation allows us to express the first two coefficients, a and b, 
in the Taylor expansion of the XeiXe2 ~^ XaXb center-of-mass frame annihilation cross 
section with respect to the relative velocity, 



cr 



xe,xe,^x^A-s t;,.ei = a + bvl, + 0{vl,) , (29) 



in terms of the partial-wave separated Wilson coefficients f^e^^<^2^^AXB^xeiXe2{^^+'^Lj^ 
and g^^i^^2-^^AXB^xe^xe2(^s+ij^^y rpj^g cocfflcient a is given by (leading order) S'-wave 
Wilson coefficients only, 

a = fCSo) + 3 fCS^). (30) 

The coefficient b receives both P-wave and next-to-next-to-leading order S'-wave Wilson 
coefficient contributions. 



( fCPi) + 3 fCPo) + fCPi) + 3 feP2 



me,+me^) V 3 3 

+ 9CSo) + l9CSo)) . (31) 



The parameters a and b in (l29l) can also be extracted numerically from computer codes 
that determine the center-of-mass frame annihilation cross sections. This is done by 
considering the cross section's behaviour for small relative velocities of the annihilating 
particle pair and performing a parabola fit to cr^^^^^^'^^^^'^ v^ei, which provides the 
corresponding coefficients a and b. Note, however, that a separation of the coefficient b 
into its constituent P-wave and next-to-next-to-leading order S'-wave contributions, as 
given in flHT]) . cannot be achieved with the sole knowledge of the cross section. Likewise, 
the separation of the S'-wave contributions for the spin singlet and triplet configurations, 
as performed in fl5U]) and flHTl) . requires intervention at the amplitude level, which is not 
straightforward for the publicly available computer codes. 

In the absence of threshold effects, resonances or enhanced radiative corrections, 
the knowledge of the coefficients a and b in X6iXe2 annihilation processes allows for 
a rather accurate calculation of the present-day relic abundance. Yet the separation 
of b into P- and S'-wave contributions is required for a consistent treatment of the 
Sommerfeld enhancement at (9(w^) because the long-range interactions responsible for 
the Sommerfeld effect depend on the quantum numbers of the incoming state. Our 
analytic approach allows us to perform this separation by construction. 
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Figure 2: Numeric comparison of the non-relativistic approximation (solid lines) to 
the tree-level annihilation cross-section times relative velocity, aWrei, for xtxt ~^ 
W~^W~^{leit) and XiXt ~^ ^ (right) reactions with the corresponding unexpanded an- 
nihilation cross section produced with MadGraph. The numeric errors on the latter 
are taken to be a Wrei/ \fN ■, where = 10^ gives the number of events used in the Mad- 
Graph calculation of each cross section value, v-^ci is given by Urei = \ve^ — Ve2 \ for the 
XeiXe2 ~^ XaXb process. The underlying MSSM spectrum is a wino-like neutralino LSP 
scenario, generated with the spectrum calculator SuSpect. The masses of the x? and Xi 
are given by m^o = 2748.92 GeV and m + = 2749.13 GeV. 

We perform a numeric check of our results for the XeiXe2 ~^ XaXb tree- level annihi- 
lation cross sections as given in fl7r|) for all initial state two-particle pairs in Tab. [T]into 
all accessible SM and Higgs two-particle final states. We consider several MSSM spectra, 
which we compute using the spectrum calculator SuSpect [23] and its implementation of 
the phenomenological MSSM, a model with 27 free parameters from soft SUSY-breaking 
terms. For each spectrum, we obtain the coefficients a and b in (!30!) and (!3T|) from our 
analytic calculation, and compare them with the corresponding coefficients extracted 
purely numerically using MadGraph [26] to calculate the cross sections. Our results 
for the coefficient a agree with the corresponding numeric expression extracted from 
MadGraph data at permille level. Similarly, we find agreement of the coefficients b de- 
rived with (13T]) and extracted from MadGraph data at 1% up to permille level, where 
the level of agreement slightly varies depending on the initial- and final-state particles. 
In addition, the level of agreement on the parameter b depends on the interval of the 
ffci variable used for the parabola fit to the MadGraph data, which for the numbers 
quoted above is taken as Vrei/c = [0,0.4]. We find that the non-relativistic approxima- 
tion is reliable for single-particle velocities up to Vejc ~ 0.3. For such velocities the 
absolute error of the non-relativistic approximation to crXeiXea-^^A^s y^^^ with respect 
to the unexpanded a-^'^i-^^a-^^-^^s y^^^ expression lies within the level of a few percent. 
Therefore the non-relativistic approximation has an acceptable accuracy for calculations 
in the early Universe during the time of Xerdecoupling, as the mean velocity of the Xci 
in that period was around Vejc ~ 0.2. 

Selected results of our numeric check with MadGraph are presented in Fig. [2] and 
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Figure 3: Numeric comparison of the non-relativist ic approximation (solid line) to avj-ei 
for the two neutral hydrogen-like two-body states X1X3 ~^ W^H~{[eit) and X1X2 
Z/"'" if "(right) to data produced with MadGraph. Again, we take the errors on the 
MadGraph data to be a Vrc\/y/N, where = 10^ gives the number of events used 
in the MadGraph calculation of each cross section value. The process on the right- 
hand side is dominated by P-wave annihilations. The underlying MSSM spectrum is 
the same as in the plots in Fig. |21 where the masses of the Xs ^ind X2 given by 
m^o = 3061.99 GeV and m^- = 3073.31 GeV. The mass of the Higgs particles takes 
the value = 167.29 GeV. 



Fig. m where the underlying SUSY spectrum contains a wino-like neutralino LSP with 
mass m^o = 2748.92 GeV and an almost mass degenerate wino-like chargino with m+ = 
2749.13 GeV. Fig. |5] shows tree- level annihilation cross sections that are relevant in the 
calculation of the neutralino LSP relic abundance including co-annihilations. The plot 
on the left-hand side displays the annihilation cross section times the relative velocity for 
the double-charged annihilation reaction XiXi ~^ W^W^ . For tVoi/c < 0.4 our analytic, 
non-relativistic approximation nicely reproduces the numeric, unexpanded cross section 
(jXi Xi ^w+w+ Furthermore, as the absolute curvature in this S'-wave dominated 

reaction is rather small compared to the coefficient a, even the absolute error that one 
would make in using the non-relativistic approximation instead of the full cross section is 
only of the order of 2% for v^c\/c ~ 0.6. The coefficient h for this reaction, calculated using 
( 13T|) . is given by 1.27 ■ 10~^^ cms. Its P- and S'-wave contributions are of the same order 
and read 2.95 ■ 10~^^cms and —1.68- 10~^^cms, respectively. The plot on the right-hand 
side in Fig. [2] depicts the single- charged annihilation reaction x^Xi ~^ tb with (massive) 
fermionic final states. As it receives significant leading order S'-wave contributions, this 
annihilation process is also relevant in the neutralino LSP relic abundance calculation 
including co-annihilation processes. Here it turns out that the b coefficient is S'-wave 
dominated, as the contributions from P-waves are suppressed by five orders of magnitude. 
Let us stress that our analytic results for the Wilson coefficients include the full mass 
dependence of the final state particles and can be applied to MSSM scenarios with fiavour 
off-diagonal sfermion generation mixing as well. 

The plots in Fig. [3] show that our results can not only be used to describe pair anni- 
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hilations of nearly mass-degenerate incoming particles XeiXe2 ~^ XaXb, but also apply 
to annihilations of a non-relativistic "hydrogen- like" XeiXe2 two-particle system of non- 
degenerate-in-mass constituents. The plot on the left hand side in Fig. |3] corresponds 
to the pair annihilation of a hydrogen-like X1X3 state into a W~^H~ final state, with 
m^o = 3061.99 GeV, which is again dominated by leading-order S'-wave contributions. 
The curvature is driven negative by the next-to-next-to-leading order S'-wave contribu- 
tions to the coefficient b, given by —5.29 ■ 10^^^ cms. The P-wave contributions are, 
however, of the same order and read 1.30- 10~^^cms. The right plot in Fig. [3]again refers 
to a hydrogen-like incoming two-body system, X1X21 where in this case the annihilation 
X1X2 ~^ H^H~ is P-wave dominated: the P-wave contribution to the coefficient b is 
given by 2.48 • 10~^^ cm s. Both the leading and next-to-next to leading order S'-wave con- 
tributions are strongly suppressed and of the order 0(10"^'^ cm^/s) and O{10~^^ cms), 
respectively. The mass of the second chargino is given by m^- = 3073.31 GeV. 

Generically, if the coefficient a in the expansion (12^ is suppressed with respect to 
the coefficient b, the curvature and hence the corresponding non-relativistic annihilation 
process is P-wave dominated. This property derives from the fact, that the leading- 
order coefficient a is related to the product of the leading order S'-wave contributions 
to the tree-level annihilation amplitude with its complex conjugate. As the next-to- 
next-to-leading order S-wave contributions to the coefficient b result from the product of 
leading order with next-to-next-to-leading order S'-wave contributions in the annihilation 
amplitudes, a suppressed coefficient a indicates a small next-to-next-to-leading order S'- 
wave contribution to the coefficient b as well. 

3.2 Analytic checks 

In [22], the authors performed a calculation of the neutralino relic abundance in minimal 
supergravity models. In the appendix, they give a complete summary of all partial-wave 
separated tree-level helicity amplitudes in XiXi ^ XaXb pair annihilations. These 
comprehensive results for tree-level neutralino LSP pair-annihilations are also referenced 
and (partly) quoted in the (SUSY) particle dark matter reviews [2] and |3], and easily 
extend to XeiXei ~^ X^Xb annihilations. Hence, these results allow for an explicit ana- 
lytic check of our expressions for the different partial-wave contributions to a neutralino 
Xe Xe ~^ XaXb annihilation cross section. The partial-wave coefficients that can be 
cross-checked in that way correspond to ^S'q-, ^Pq-, ^Pi- and ^P2-wave XeiXei ~^ XaXb 
annihilation reactions, and the leading order and next-to-next-to- leading order ^S'o-wave 
contributions can be compared separately. As already inferred from ([8]) and noted at 
the end of Sec. 12. 2[ there are no ^S'l and ^Pi partial- wave contributions for annihila- 
tion reactions of identical incoming particles, which is the case covered by [22]. Our 
expressions for the partial- wave separated XeiXei ~^ XaXb annihilation cross sections 
into all possible SM and Higgs final states, obtained from fl27|) . agree with the corre- 
sponding terms derived from the helicity amplitudes in \2!2y\ . We note that our results 

^ The only minor discrepancies that we find are related to P-wave contributions: our results for 
■^Pi-wave XeiXei ~^ H~^H~ annihilations correspond to a factor 2 instead of a factor 4 in the second 
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for annihilations into a pair of fermions include the case of flavour-off-diagonal sfermion 
mixing as well, which is covered in [2] and [3], but was not yet included in [22], wherein 
only flavour-diagonal right-left sfermion mixing was taken into account, although it is 
straightforward to extend these results to the general flavour-off-diagonal case. 

The comparison with analytic results for inclusive leading-order ^5*0- and ^S'l-wave 
pair-annihilation reactions of a pure wino-like neutralino Xi ^^^^ its niass degenerate 
chargino partners xf i^ito all possible SM and Higgs final states considered in Ref. [H] 
provides another useful check of our results for the absorptive part of the Wilson coef- 
ficients. The results in [13] comprise all possible neutral, as well as single and double 
charged inclusive pair-annihilation reactions. The masses of the SM and Higgs particle 
final states are set to zero, such that the corresponding results can be understood as the 
leading-order term in an expansion in msM/iTiy^o and mmggs/'mxl- Furthermore, all super- 
symmetric particle states heavier than Xi and xf treated as completely decoupled. 
We agree with all results for the inclusive annihilation reactions given in [2]. In partic- 
ular we agree with the results in [T3] that refer to leading-order ^S'o-wave XiXi ~^ XiXi 
as well as XiXi ~^ XiXi reactions, which can be related to the Wilson coefficients 
fXiXi^Xi Xi (^Sq) and /^i So) , therewith permitting an explicit check of some of 

our Wilson coefficients encoding off-diagonal scattering reactionso 



4 Discussion 

4.1 Unitary vs Feynman gauge 

The computation of the absorptive parts of the Wilson coefficients for forward-scattering 
reactions, XeiXe2 ~^ ^a^b ~^ XeiXe2J has been performed using both the unitary and 
Feynman gauge. The results agree numerically, which provides a further check of our 
calculation. For the off-diagonal reactions, where the incoming and outgoing states are 
different, the use of unitary gauge for final states with two massive vector bosons in 
the final state introduces enhanced 1/My and 1/My terms proportional to the mass 
differences between the incoming and outgoing particle species, which must cancel in 
the final result. Similarly, a cancellation of 1/My enhanced terms in off-diagonal rates 
with one massive vector boson in the final state has to take place. However, for these 
cancellations to occur, one has to also expand the SUSY mixing matrices systematically 
in the gauge boson masses My- In the same way, the mass differences between the 

term of Eq. (A27b) in [22]. In the case of '^Po-wa.ve XeiXei ~^ ff reactions, our results correspond to a 
factor ■\/2/3 instead of a factor \/6 in the second term in the first hne of Eq. ( A29b) in [22] . 

^ The authors of Ref. [TJ] also provided analytic results for exclusive leading-order -'^S'o-wave anni- 
hilation reactions for both the cases of a wino-like and a Higgsino-like neutralino LSP scenario in a 
previous work |13| . We agree with the results for all diagonal XeiXe2 ~^ XeiXe-i reactions. A typo in the 
off-diagonal terms in Eq. (28) of [T3] was fixed in [T3|, and the latter agrees with our findings. In the 
Higgsino-like scenario, we get differing expressions for off-diagonal XiX? ~^ XxxX xTXi^ ~^ XiX? 
reactions in the case of W^W~ and ZZ final states for both i = 1,2: our results are a factor 4 and a 
factor 2 larger, respectively, than the corresponding expressions presented in [13] . 
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incoming and outgoing particles have to be expanded in My and in the differences of 
soft SUSY breaking parameters Mi, M2, /x, if these differences are small. The latter 
expansions must be done differently depending on how many neutralinos and charginos 
are (nearly) mass-degenerate. The presentation of the results computed with unitary 
gauge then has to distinguish among many cases and also consider diagonal and off- 
diagonal terms separately, since for the diagonal terms it is desirable to keep the full mass 
dependence as well as unexpanded mixing matrices. We thus find it more convenient 
to use Feynman gauge for the calculation of the off-diagonal reactions, which allows to 
keep the coupling matrices unexpanded and a more concise presentation of the results. 
The price for this is that one must compute a large number of unphysical final states 
containing pseudo-Goldstone Higgs and ghost particles, see Tab. [3l 

4.2 Off-diagonal terms 

Our framework aims to describe the annihilation of a pair of non-relativistic charginos 
or neutrahnos {xiXj) iiito SM and light Higgs particles pairs [XaXb) including potential 
interactions between all nearly mass-degenerate XX states, that can produce a Sommer- 
feld enhancement of the rates. A contribution to these enhanced annihilation rates is 
given by the imaginary part of the amplitude for a process of the type, 

XiXj • • • XeiXea ^aXb ^ Xe^Xez ^ ... ^ XiXj , (32) 

where the intermediate states involved in the short- distance annihilation, XeiXe2 ^ind 
Xe4,Xe3, can be different (off-diagonal annihilation terms), compare to Fig. [T]for a figu- 
rative illustration. In a recent work [TS], a general formalism which also aims to com- 
pute the Sommerfeld-enhanced annihilation rates for a coupled system of neutralino and 
chargino pairs, has been presented which, however, does not implement the possibility 
of off-diagonal transitions in the hard part of the annihilation process. We show in this 
section that the off-diagonal terms can indeed be relevant, and should be accounted for 
in the calculation of the Sommerfeld enhanced rates. 

Naively, if the final state Xa_Xb is allowed for both XeiXe2 and XetXes, given one par- 
ticular partial-wave configuration of the two-body systems, the off-diagonal absorptive 
amplitude can be of the same size as the diagonal absorptive amplitude, i.e. 

JldF^As] AiXe^Xe2 XaXb) X AiXe.Xe, ^ XaXb)* 

~ j [dPS ab] \A{Xeae2 ^ XaXb)\' , (33) 

since the phase-space integration involves very similar kinematics. An example of the 
latter is given by the annihilation rates of wino-like neutralino dark matter, where the 
wino-like neutralino (x?) is highly degenerate with its charged SU{2)l partners (xf )• 
In such scenario the spin-0 XiXi system mixes with the XiXt state through VT-boson 
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exchange. The inclusive annihilation rates that have to be fed into the calculation of the 
enhanced rates for the spin-0 XiXi channel in the wino limit read 



^ /{11}{11} I '^oj 



(34) 
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where 0^2 = S'l/^^? 92 denotes the SU{2)l gauge coupling, and all gauge boson and 
Higgs-particle masses are treated as massless. We see explicitly that the off-diagonal 
term (!36|) is of the same order as the diagonal reactions (134] - [35|) . 

In order to stress the importance of the off-diagonal annihilation terms, we have 
computed the thermally averaged effective annihilation cross section, (cTcfrf), which en- 
ters the Boltzmann equation for the calculation of the dark matter yield, for the same 
wino-like scenario used for the checks with MadGraph presented in Sec. 13. ![ and com- 
pared to the results obtained when the off-diagonal terms are switched off by hand. 
For the necessary formulae to compute ((Tefft'), including co-annihilation effects, we re- 
fer the reader to j27l[2S] . The annihilation rates have been calculated using fl27|) with 
Wilson coefficients multiplied by the appropriate Sommerfeld factors computed solving 
the coupled-channel Schrodinger equation for each partial-wave. The details about the 
calculation of the Sommerfeld enhancement factors from the long-range interactions will 
be given in a future publication p¥j . 

We observe from Fig. H] that removing the off-diagonal annihilation terms decreases 
the thermally averaged cross section by a factor larger than 1.5 at small temperatures. 
The corresponding thermal relic abundance of the dark matter in the present Universe, 
^'Duh'^i obtained by numerical integration of the Boltzmann equation, gets then increased 
by approximately 20%, if the off-diagonal reactions are neglected. The latter represents 
thus a sizeable effect which has to be accounted for in such a scenario. 

5 Summary 

The calculation of the thermal relic abundance of the lightest neutralino as a promising 
dark matter candidate within the MSSM places strong bounds on the MSSM parameter 
space, assuming that the observed cosmic dark matter has particle nature and is com- 
posed solely of the neutralino LSP. Given the expected future experimental accuracy of 
the measurement of the cosmic dark matter abundance observed today, radiative cor- 
rections to the pure tree-level annihilation cross section, entering the relic abundance 
calculation as a central ingredient, should eventually be taken into account. The in- 
clusion of 1-loop corrections to the annihilation cross section as well as the systematic 
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x = m. ^IT 

Figure 4: Thermally averaged cross section as a function of x = vti^jT with fixed = 
2748.92 GeV, including the Sommerfeld effect (solid blue line). The same quantity 
computed with the off-diagonal perturbative annihilation rates set to zero is depicted by 
the dashed red line. The perturbative result is also shown as a dotted line for comparison. 
The decrease of (cJefrw) towards a; ~ 1 is due to large negative C(p^) terms in the XX 
annihilation rates, which are unphysical because the non-relativistic expansion (p7|) for 
the annihilation rate becomes unreliable for large temperatures. 



treatment of Sommerfeld enhancements, has recently been a field of elaborate studies 
in the literature. Similarly, in the context of dark matter annihilation processes in the 
present Universe relevant in indirect detection, the above types of radiative corrections 
to the neutralino pair-annihilation cross section have been studied extensively. 

In this paper we take advantage of the non-relativistic nature of the annihilating 
neutralinos in the present Universe as well as during thermal dark matter decoupling 
in the early Universe, which introduces a clear separation of energy scales in all anni- 
hilation processes of interest. The latter property allows us to set up an effective field 
theory (the NRMSSM) of non-relativistic neutralinos and charginos, that provides an 
appropriate setup for a systematic investigation of radiative corrections to neutralino 
LSP pair annihilation processes both in the present and the early Universe, taking co- 
annihilations with nearly mass degenerate neutralinos and charginos into account. As 
a first step in the explicit construction of the NRMSSM we have derived fully analytic 
formulae for the absorptive part of the Wilson coefficients of four-fermion operators in 
the effective theory pertaining to S-wave annihilation, that encode the hard annihilation 
rates of XeiXe2 ~^ XaXb processes (see ([6]) and Tab. [2]). Our results separately include 
leading-order ^5*0- and ^Si-wave as well as all P-wave and next-to- next-to-leading order 
S'-wave Wilson coefficients and apply to general neutralino and chargino states in the 
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MSSM. Flavour off-diagonal sfermion generation mixing can be covered, and we k;eep the 
full mass dependence of all SM and Higgs particles. Analytic results for the absorptive 
part of leading order S-wave Wilson coefficients are presented in the appendix. Results 
for P-wave and next-to-next-to-leading order S-wave coefficients will be given in a fu- 
ture publication [23]. By taking into account charge-neutral annihilation processes of 
a chargino pair as well as singly charged and doubly charged annihilation reactions of 
non-relativistic neutralinos and charginos, we extend the analytic results for partial wave 
decomposed neutralino LSP pair-annihilation cross sections given in the literature |22] . 

We have shown that the non-relativistic expansion to produces accurate results 

up to Vrei ~ 0.6, which is sufficient for relic density computations, and certainly for 
dark matter annihilation in the present Universe. Our analytic results may therefore 
substitute for time-consuming numerical computations. 

Our aim is to apply the effective field theory formalism to the calculation of Sommer- 
feld-enhanced (co-) annihilation cross sections in the neutralino relic abundance calcu- 
lation. As scattering prior to the annihilation process can lead to transitions from an 
incoming particle pair to another nearly mass degenerate neutralino or chargino two- 
particle state, a proper treatment of the Sommerfeld effect requires the knowledge of the 
absorptive part of off-diagonal annihilation rates, XeiXe2 XaXb — > XeiXe-s^ for ah pos- 
sible SM and Higgs two-particle states XaXb (see Fig. [1]). To the best of our knowledge 
we present for the first time analytic results that allow for a systematic treatment of all 
these off-diagonal rates in Sommerfeld-enhanced (co-) annihilation reactions for general 
masses and composition of the Xci particles. The implications of these results for MSSM 
relic density calculations will be studied in a forth-coming publication |24j . 
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A Absorptive parts of Wilson coefficients of dim- 
ension-6 operators in (5£ann 

We present the leading order contributions to the absorptive part, jxx^xxj^ss+i^^^^^ ^j^g 
Wilson coefficients that correspond to the local four-fermion operators given in Tab. |2l 
The /^'=i^'=2~^><:<i4><:<!3 (^2s+i^^^ encode the absorptive part of hard 2 — 2 scattering reactions 
of an incoming particle pair XeiXe2 of non-relativistic charginos or neutralinos in a given 
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XX XX 


vs ss ff r]f] 


X'X^ 

A A A A 


^^'"^^ G^G- ,G^H- , v^v^ ^ ^ 
W-G+,W-H+ H+G-,H+H-, 




ZG+,-iG+, G+h°,G+H^, 7]+f]^, 
W+Z, ZH+,jH+, H+h°,H+H^, u^d^, r]^f]-, 
W+-f W+h^,W+H°, G+G'^,G+A^, z/V r/+r/^, 
W+A^ H+G°, H+A^ r]^f]- 




^ ^ G+G+, 
W+G+ 

W+W+ ^ ' G+H+, 
W+H+ ^ ' 
H+H+ 



Table 3: Particle pairs XaXb in XX ~^ XaXb — ?• XX scattering reactions (abbrevi- 
ated as XX ~^ XX)-, that we account for in the calculation of the absorptive part of the 
Wilson coefficients, classified according to their type: VV, VS, SS, ff and rjfj. Nega- 
tively charged processes, corresponding to the charge-conjugates of the singly or doubly 
positively charged reactions above are not explicitly written. 



'^^+^Lj partial- wave state into an outgoing non-relativistic Xe4Xe3-pah in the same partial- 
wave configuration. They allow to reproduce the inclusive tree-level center-of-mass frame 
annihilation cross sections of a non-relativistic XiXj'PaiJ^ i^ito SM and light Higgs two- 
body final states XaXb, expanded in the relative velocity of the annihilating particle 
pair. The general case includes off-diagonal processes XeiXe2 ~^ XaXb — )■ Xe^Xea with 
XeiXe2 7^ XeiXes, foT all pairs of non-relativistic neutralinos and charginos. Since the 
jxx^xx^^s+^l^j'j gj-Q infrared-safe at leading order, we are able to give analytic results for 
the individual contributions f^^~^^A^B-^^^{^^^'^Lj) pertaining to an exclusive final state 
XaXb- 

A.l Notation and definitions 

Recall that the calculation is performed in Feynman gauge. Hence the two-particle final 
states XaXb that we account for can be classified to be of vector- vector (V^^), vector- 
scalar (^5*), scalar-scalar {SS), fermion-antifermion (//) or ghost-anti-ghost {rjff) type. 
They are listed in Tab. [3l The determination of the absorptive part of the Wilson 
coefficients for the processes XeiXe-i XaXb Xe^Xe^ requires the calculation of a 

■^The covered XiXj-states have been collected in Tab. [H 
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selfenergy (ss) 



Figure 5: Generic selfenergy- diagram in XX ^ ^a^b XX reactions. Particles Xa 
and Xb represent any two-body final state of SM and Higgs particles, which can be 
produced on-shell in xX ~^ XaXb annihilations. 



large number of Feynman diagrams. To be able to present the results in an efficient 
manner it is convenient to make use of the classification in VV-, VS-, SS-, ff- and rjf]- 
type XaXb particle states and to further subdivide the contributing diagrams according 
to their topology. In each of the classes under consideration there arise generic 1-loop 
amplitudes with selfenergy-, triangle- and box-topology shown in Figs. [SHHl The 
generic selfenergy-diagram as well as the four generic triangle- and box-diagrams cover 
all possible kinematic configuration^ that can arise in a XeiXe2 ~^ XaXb — >■ Xe4Xez 
1-loop amplitude. Note that we have assigned specific directions for the fermion flow in 
each diagram in Figs. [SHHl indicated by the arrows, as it is convenient in the context 
of calculations involving both Dirac and Majorana fermions, following the Feynman 
rules for fermion-number violating interactions set out in [32]. The depicted fermion 
flows establish our convention to arrange the external fermion states Xei! ^ = • • • > 4 in 
descending order, see Tab. |2i 

We calculate analytically the absorptive part of any of the contributing selfenergy-, 
triangle- and box- amplitudes, subject to our convention for the fermion flows. Thereby 
we consider generic external Majorana fermions, generic t- and u-channel exchanged Ma- 

''The case of four different triangle- and four different box-diagrams in Fig. [6] and Fig. [7] applies to 
non-identical particles Xa Xb- For identical particles Xa = Xb, triangle (box) 1 and 3 as well as 
triangle (box) 2 and 4 coincide. In this case only one of the identical diagrams must be taken into account 
to compute the corresponding Jxx^XaXa^xx coefficients. This rule incorporates the symmetry factor 
of 1/2 for identical final-state particles, that one would take into account in the conventional calculation 
of the tree-level XeiXe2 ~^ XaXa annihilation rate. 





triangle 1 {tis) 



triangle 2 (sti) 



triangle 3 {t2s) 



triangle 4 (5^2) 



Figure 6: Generic triangle-diagrams in xX ~^ XaXb XX reactions. 
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box 1 (tita) box 2 [titi) box 3 (^2^1) box 4 (^2^2) 



Figure 7: Generic box-diagrams in xx ~^ XaXb — > XX reactions. 



jorana fermions or sfermions, generic XaXb states of type VV, VS, SS, ff and rjr], and 
hence use generic 'place-holder' coupling factors at each vertex. This allows us to deter- 
mine the generic form of those terms in the contributions to the f^^^^A^B^xx^'^s+ij^^^^ 
that are associated with the kinematics of the xX ~^ XaXb — )■ XX reaction, where each 
of these kinematic terms multiplies a certain combination of the place-holder coupling 
factors. In particular, these kinematic contributions are generic in the sense that they 
apply to both the cases of external and internal Majorana and Dirac fermions. 

A specific diagram's contribution to the absorptive part of a particular XeiXe2 ~^ 
XaXb — XeiXe^ MSSM 1-loop process is obtained by replacing the generic place- 
holder coupling factors with their actual expressions in the above described generic 
Majorana fermion 2 — ?■ 2 scattering reactions. Note that by choosing these coupling 
factors properly, all XeiXe2 ~^ XaXb — )■ Xe^Xez processes with external and internal 
Majorana or Dirac fermions can be covered, although the kinematic contributions are 
calculated referring to the generic Majorana fermion 2 — >■ 2 scattering reaction. Hence, 
the absorptive part of the Wilson coefficient, which encodes the absorptive part of a 
XeiXe2 ~^ XaXb — )■ XeiXe^ scattering reaction, with the incoming and outgoing two- 
particle states in a "^^^^Lj partial-wave configuration, can be written as 

{eie2}{e4,ez} V J) 



Xb V Xa Xa 




box 1 {tit2) box 2 {titi) box 3 (^2^1) box 4 (^2^2) 

Figure 8: Generic box-diagrams in xX ~^ XaXb XX reactions, with XaXb a pair of 
SM fermions. 
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\ ri 11,12 
4 

, \^ (Q!)XeiXe2^^A^i3-)>Xe4Xe3 ^{ci) X aX b /2s + l T 

+ Z^Z^Z^'^n,ni2 *-^n,ni2 I 

Q!=l n 41,12 



= 1 n 21,42 



Here 0^2 = qH^t^^ where §2 denotes the SU{2)l gauge couphng. The sums in the first 
hne on the right-hand side of fl371) collect all contributions from selfenergy-amplitudes. 
Similarly, the second (third) line gives the triangle- (box-) amplitudes' contributions. We 
use the index a to enumerate expressions related to the four different triangle- and box- 
amplitudes0 according to the labeling of the diagrams in Figs. EHSi Further, we indicate 
the kinematic factors of the generic 2 — 2 Majorana fermion scattering amplitudes 
within a given class and topology with capital letters (i?„ j^j^i ^n"iii2' -^n"iii2)- These are 
the quantities that include the kinematics of the process and hence encode the ^^~^^Lj 
partial-wave specific information. The process-specific coupling factors that multiply the 
kinematic factors are denoted with lowercase letters (^„,iii2' '^n"!i42' '^n"iii2)- Depending 
on the type of the particles Xa and Xb as well as the topology, there is a fixed number 
of different coupling-factor expressions that can occur, together with the corresponding 
kinematic factors. The different contributions are enumerated with the index n in fl37|) 
above. Finally, in each of the processes there is a certain set of particle species that 
can be exchanged in the s- or the t-channels of the contributing amplitudes. These are 
labeled with the indices ii and 22- 

The generic structure of the Wilson coefficients in ( 137|1 suggests to give the coupling 
factors and the kinematic factors separately. A recipe for the construction of the coupling 
factors jjjj, c^"lii2' '^n"iii2 covered reactions is given in Sec. IA.2I Analytic 

results for the kinematic factors j^^^, C^l_^^^, d'^^\_^^^ for the leading-order ^5*0 and ^5*1 
partial-wave configurations can be found in Sec. IA.3I These expressions depend on the 
masses of the external and internal particles in a particular XeiXe2 ~^ ^a^b Xe^Xez 
process. However, the kinematic factors are generic in the sense that their form is 
the same for all possible external two-body states XeiXe2 and XezXen of neutralinos or 
charginos and all XaXb particles within one of the classes W, VS, SS, ff or rjf]. 

The coupling and kinematic factors will depend on the supersymmetric particles' 
mixing matrices and masses, respectively. We adopt the same notation as in [33] and 
hence introduce the chargino and neutralino mixing matrices Z± and Z^r defined via 



Z+ = I ""^^ I , (38) 
^x^ 



^For identical particles Xa = Xb the index a has to be taken from 1 to 2 only, see footnote ID 
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ZnM^o 



V 



A2 



(39) 



where M^± and M^o denote the chargino and neutrahno mass matrices, respectively (for 
details regarding the mass matrix expressions refer to [33]). m^+,j = 1,2, and m^o,i = 
1, . . .4 indicate the masses in the mass eigenstate basis of charginos and neutralinos. 

In order to properly apply the formulas for coupling and kinematic factors in Sec. IA.2I 
and Sec. I A. 3 1 given a specific MSSM spectrum, it is important to note that the NRMSSM 
and hence the analytic expressions for the Wilson coefficients explicitly rely on the posi- 
tivity of all mass parameters. This derives from the fact that the NRMSSM Lagrangian 
is obtained by extracting the high-energy fluctuations (of the order of the particle mass) 
from the relativistic fields, which yields the non-relativistic kinetic term £kin shown in 
(|3]). For species other than the LSP, the procedure leads to the mass- difference terms 
{rrii — mLSp) in ([3]). If any of the rrii in £kin. is negative, then the corresponding mass 
difference counts as 0(mLSp), an indication that the parameterization used to relate the 
relativistic and non-relativistic fields for that particle species is not the appropriate one. 
The simplest way to obtain the NRMSSM Lagrangian in case that the mass m^^ of 
one or several of the external Xei particles happens to be negative for a given MSSM 
spectrum, is to perform a field redefinition of the corresponding MSSM fields that yields 
mass terms with positive mass parameters. Such a field redefinition affects the chargino 
and neutralino mixing matrices, which are mapped in the following way: 



( 



'sgn(m^4 



'sgn(m^,o^ 




sgn(m 



sgn(m^o^ 



(40) 



sgn(m^o, 



sgn(m^o 



(41) 



xV ) 



(We define -y/— 1 = i.) The redefined mixing matrices Z± and as well as the corre- 
sponding positive mass parameters for all MSSM neutralino and chargino fields should 
be used within the expressions given in Sec. IA.2I and Sec. IA.3I 



A. 2 Coupling factors 

By construction, the absorptive part jx^\X^i^^AXB^xe^x^-i, Qf individual Wilson co- 



efficient is associated with the product A 



(0) 

Xe^X^^^^A^B 



X {A. 



(0) 



)* of Born- level 
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diagram ti diagram t2 diagram s 

Figure 9: Generic tree- level amplitudes in XX ~^ X^Xb annihilations, referring to 
VV, VS and SS-type final state particles X^Xb- The generic form of s-channel exchange 
diagrams for X^Xb = Tjfj final states agrees with the s-channel diagram above. The 
vertex-factors V^^^'^^ are defined as Vj^^'^'' = "J^ir^^^ + Q'ef 75), if attached to a three-point 
vertex with a gauge-boson (with Lorentz-index p), and Vf^'^^ = {rfi +qei'^^)i if associated 
with a vertex that involves a scalar particle X^, Xb or X^. Here the expression Tf^{q^^^) 
either denotes a vector or scalar (an axial- vector or pseudo-scalar) type of coupling factor. 
For the definition of CABXi and the Lorentz structures Labx^ we refer to Tab. S] below. 



annihilation amplitudes A^^^ related to XciXej XaXb reactions, integrated over the 

XaXb two-particle phase space, see (fT6l) . Each of the tree-amplitudes -^xx-^XaXb 
ceives contributions from diagrams with t-channel neutralino or chargino exchange as 
well as from diagrams with s-channel Higgs-particle or gauge-boson exchange, such as 
the generic diagrams shown in Fig. [91 In case of fermionic final states XaXb, instead 
of neutralino or chargino t-channel exchange, t-channel sfermion-exchange occurs, as 
depicted in Fig. [TUJ Note, that in Fig. M and Fig. [TU] we again have established a spe- 
cific fermion flow, which in particular coincides with the convention for the fermion flow 
associated with the incoming two particles in the 1-loop amplitudes in Figs. [SHU 

A contribution to the amplitude A^^^_^XaXb iii^olves a product of two coupling factors, 
coming from the two vertices in the tree-level diagrams. The generic form of these vertices 
is indicated in Fig. |9] and Fig. [TOl It is especially convenient to write all vertex factors in 
any of the amplitudes contributing to the non-relativistic xx ~^ XaXb — )■ XX scattering- 
processes as a combination of (axial-) vector or (pseudo-) scalar coupling factors, instead 
of using left- and right-handed couplings, as it is common in calculations related to the 
MSSM. The reason for that is, that in the non-relativistic limit, either the contributions 
to the annihilation amplitudes involving the axial-vector (pseudo-scalar) coupling will be 
suppressed with respect to the corresponding contributions related to the vector (scalar) 
coupling, or vice versa, such that the use of (axial-) vector and (pseudo-) scalar couplings 
allows for a clearer understanding of leading and suppressed contributions in the non- 
relativistic scattering regime that we aim to study. 

Each of the coupling factors 6„, ci"'* and rfi"^ that occur in (EZ]) is given by a product 
of two coupling factors, r or g, arising in an individual diagram in ^XaXb^ ^^"-^ 
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diagram ti diagram t2 diagram s 



Figure 10: Generic tree-level amplitudes in xx ~^ X^Xb annihilations, with X^Xb = 
f^f^. For the definition of V.f '^^ see Fig. O The generic vertex factor V^^l^ is defined 

Vxli = ''^i^L + Q^Keil^i s^ch ^h^^ ''^xli i^^Kei) denote coupling factors of scalar 
(pseudo-scalar) type. 



the complex conjugate of another such two-coupling factor product originating from 
.^XaXb' "'■^ following, we give a recipe how to construct the coupling factors 
in ( 1571) for a specific XeiXe2 ~^ XaXb — )■ XciXea reaction, such that taken together with 
the kinematic factors in Sec. IA.3t they allow to determine the absorptive part of the 
Wilson coefficients /: 

1. Draw all tree-level diagrams that contribute to XeiXe2 ~^ XaXb and Xe^Xe-j, 
XaXb annihilation amplitudes, analogous to the generic diagrams sketched in 
Fig. [9] or Fig. [TOl In particular, assign the same fermion fiow as indicated for the 
generic diagrams. 

2. Determine the process-specific (axial-) vector and/or (pseudo-) scalar coupling fac- 
tors, that arise instead of the generic q^^f or r^^f place-holder expressions at the 
generic amplitudes' vertex factors. As the xX ~^ X^Xb processes may involve Ma- 
jorana as well as Dirac fermions, and the latter involve a conserved fermion-number 
fiow, note the following rules: 

a) If the direction of the fermion-number fiow related to a Dirac particle coincides 
with the direction of the fermion fiow (fixed as in Fig.[9]and Fig. fTOl) . the xX ~^ 
XaXb process specific coupling factors at the vertices are directly deduced 
from the corresponding interaction terms in the underlying Lagrangian. These 
coupling factors are given later in (144] - [50|1 . 

b) Otherwise, if the fermion-number fiow is antiparallel to the indicated fermion 
fiow, vector coupling factors at vertices attached to a Dirac fermion line, are 
given by a factor —1 times the expression for the vector coupling given in (jUJ- 

Axial- vector, scalar and pseudo-scalar coupling factors are unchanged 
with respect to case a) above. 
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3. Build all possible two-coupling factor products, including possible signs related to 
vector couplings, as far as the case in 2b) above applies, that can arise in each 
single diagram. 

4. Multiply each of the two-coupling factor products, that arise in the A^^^ ^XaXb 
amplitude, with the complex conjugate of each of the two-coupling factor products, 
arising in A^^^^ ^XaXb' ^ result, all possible coupling factor combinations that 
can occur in f ^<^i^<^2^^AXB-^xe^xe3 gj-Q obtained. 

Rule 2b) arises in the following way for the case of diagram s in Figs. [QlfTOl according 
to our convention for the fermion flow in Fig. [9], we obtain an expression —v{pi)ru{p2) 
for the incoming particles' spinor chain if the case under 26) applies, where F denotes 
the involved Dirac-matrix structure. The minus sign accounts for our convention for the 
order of the external fermion states. This expression can be rewritten as 

-v{p,)Tu{p2) = v{p2) C F^C-^(pi) , (42) 
wherein C denotes the charge conjugation matrix. Using 

C ^ I -r for F = 7. , ^^^^ 

[ F for F = 1,75,7^75 , 

the origin of the minus sign rules for vector couplings under 2b) above becomes obvious. 
For diagrams with t-channel exchange, a similar derivation also confirms rule 2b). 

Let us introduce the shorthand ad to indicate the diagrams a and d in the XeiXe2 ~^ 
XaXb and XciXe-i ^a^b processes, respectively, to which the coupling factors in a 
specific coupling factor combination are related. Both a and d can be given by s, ti or t2, 
see Figs. [914T01 Coupling factor combinations originating from ss lead to the b factors, 
that correspond to the generic selfenergy-amplitude in Fig. O We label coupling factor 
combinations, that originate from tis, sti, t2S and with the superscript a = 1, ... ,4, 
respectively. These coupling factor combinations, related to one t-channel and one s- 
channel exchange diagram give rise to the c^"-* expressions in fl37|) . The a = 1 , . . . , 4 
label-convention for the specific coupling factor combinations allows to correctly allocate 
the c^"^ to their corresponding generic triangle-amplitude 'triangle a' in Fig. Ei Coupling 
factor combinations originating from tit2, ti^i, t2ti and t2t2 are labeled with superscript 
a = 1, . . . , 4, and give rise to the ci^"-* expressions. As in case of the c^"\ this convention 
correctly assigns d^""^ expressions to their corresponding 'box a' amplitude in Fig. [7] or 

Fig. [a 

We introduce the index n in order to label the different coupling factor combinations 
for a given fixed ad. Each individual n is given by a character-string, where the ith 
character gives the type (r or q) of the coupling factor which is related to the ith vertex 
in the particular diagram ad in Figs.|5HHl Thereby we have enumerated the vertices in the 
respective 1-loop amplitudes from top to bottom and left to right. Vertex factors of type 
CABXi are not specified in the corresponding string n, because the nature of the particles 
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ig2 b'l i^j + QtjTo) 



igi [7'1 {r^j + qijlb) 



192 [Shj+Plijlb) 



Figure 11: Generic form of the vertex factors in three-point interactions of neutrahnos 
and charginos with SM and Higgs particles, upon which our definition of the (axial-) 
vector and (pseudo-) scalar coupling factors given in the text is based. 



Xai Xb and Xi involved in the diagram completely characterize that coupling. For 
triangles with X^Xb = VV, VS or SS, for example, the index n will range over strings 
rrr,qqr, where the characters r or g indicate the type of coupling of the external 
particles to the X^Xb pair and to the single s-channel exchanged particle species, see 
Fig. El 

To further specify the coupling factor combinations for given aa and n, we use the 
labels ii and 22 to indicate the particle species that are internally exchanged in diagrams 
a and a. Therewith, the coupling factor combinations fcniii2! c^^iiia ^n'iii2 ^^^^ should 
enter in ( 137]) . together with the generic kinematic factor expressions given in App. lA.St 
can be unambiguously determined. 

In order to completely fix our conventions, we summarize in the following the ex- 
pressions for the (axial-) vector and (pseudo-) scalar coupling factors, that arise in the 
three-point interactions of charginos and neutralinos with SM and Higgs particles. The 
definitions of the couphng factors assume that we take xt to be the particle and x7 its 
anti-particle, such that the Dirac fermion number flow, indicated by the arrow on the 
Dirac fermion line for a chargino, will always refer to the direction of xt flow. The latter 
convention agrees with that of Rosiek |33] • 

The generic form for the vertex factor, that describes the 3-point interaction of an 
incoming neutralino Xp outgoing chargino xt either an incoming charged Higgs 
particle or or an incoming """-boson is given in the left-most diagram in Fig. [TTJ 
Note that the gamma matrix 7^ in the vertex factor has to be omitted in case of inter- 
actions with the charged Higgs particles. The specific scalar and pseudo-scalar or vector 
and axial-vector coupling factors, that have to be replaced for the generic r^j and qij 
couplings therein read 



{p. 



w 

V 



( Zfj*Zl'* + —Zl'*{Z^J* + ta.ndwZ 



N 



± Z]^ [Z^j^Z^: - ^Z^-iZ'J + tan^H/4'; 



2i* 

+ I ' 
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7! ^ -^7! ^ + 



(44) 



where i//' = and ifg*^ = , and the mixing matrices are defined as in Ref. [33] . 
The generic form of the three point interaction of either two neutrahnos or two charginos 
with an electrically neutral gauge boson or Higgs particle is depicted in the second dia- 
gram in Fig. [TT] Again, the gamma-matrix 7^ has to be omitted in the vertex factor if 
the corresponding reaction refers to interactions with the neutral Higgs particles. In the 
case of an incoming Xj and an outgoing xt ^ ^^e (axial-) vector and (pseudo-) scalar cou- 
plings, that encode interactions with the neutral scalar and pseudo-scalar Higgs particles 
[h^, H^, G^, A°), the Z-boson or the photon are given by the following expressions: 



rrO CtO 



1 



4r (pfr 



2V2 

i 

2V2 



^Im ^ ^2j* ^li* 



-^(ZHZ'J* + Zl'*Zl' + 2{cl 

4:CW V 



W 



4ciy 



li ylj* 



Z''Z 



(45) 



where = H'^, if 2 = ^^'^ = A^, A2 = G^. Finally, three-point interactions of an 
incoming and an outgoing Xi with a (pseudo-) scalar Higgs particle or the Z-boson 
involve the following (axial-) vector or (pseudo-) scalar coupling factors 



r72m n'ii* ryXm ry'M* 



ry2m nrii* rylm ryZi* 



Z'^ — tan 6w Z^fj 



w 



Zfj* — tan 6w Zl^j* ] + (i ^ 



w 



3} 



4c 



w 



The (axial-) vector and (pseudo-) scalar coupling factors in fHS]) and 
related to interactions with neutral SM and Higgs particles, satisfy 



± c.c. , 
± c.c. , 
(46) 

which are all 



V. 



V 



* * 

^ij ~ ^ji 5 Pij ~ ~Pji ■ 

as a consequence of the hermiticity of the underlying SUSY Lagrangian. 



a 4 



(47) 
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The generic form of the vertex factor for three-point interactions of a neutrahno or 
chargino with a SM fermion and a sfermion is given in the right-most diagram in Fig. [TTJ 
In case of interactions of an incoming SM fermion with a sfermion 5"^ and an outgoing 
neutrahno Xi? the specific (pseudo-) scalar couphngs, that have to be replaced for the 
generic Snj and puj expressions in Fig. [H] read 

db(^dD\ _ 1 ^ 7l»7(/+3)j ( y3» y/j I y-ii v(I+3)j 

Iij[Plij) - Qd lanUw^N " T—E: TT:;— \ ^ ^N^D 

V2 2v2cos/3Mvi/ ^ 



s 



T ^ (t,Z^ + (g, - Trf)Z^^ tan^H/) 4' , 
= ^ ^ {^'^ - Z'n^^^(^w) Z'r , 

s%{p%) =^q, tan^vK^il"4'^'^' - { Zf^*Zi' ± ^^^4^+^^^' 

Rjyt^ivJ ^2 ^' TV L 2v^cos/3MvkV ^ ^ ^ ^ 

T (r^Z^* + {qi - Ti)Z]i tan Z[' . (48) 

/ = 1,2,3 denotes the generation index for the fermions, and j = 1, ... ,6 labels the 
sfermion states {j = 1, 2, 3 in case of sneutrinos uj). Tj and qj are defined as 

Tu = —Td = —Ti = - , 

2 1 , , 

Qu = - , qd = -- , qi = -I . (49) 

The superscripts, fS, on the couplings in f HS]) . refer to the fermion (/)- and sfermion (S)- 
type involved in the underlying interaction. In case of chargino-fermion-sfermion inter- 
actions we find (a sum over repeated indices is implicit) 

T \k''*Z'_^Z'J , 

sf^lpfJ) = ^ K''z''*zi^* ± K'^Z'^zif^'^'* 

Rjyi^irjJ 2V2 cos (3 Mw " 2V2 sin (3 Mw 
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LABXi 




y \rvA i^B ) y \i^t i^A) 


^ y K'^B f^'i ) 








v^ss 


{ks - KY 




ssv^ 


{ks - kAY 




sss 


Mw 






^B 




vv s 


Mw 





Table 4: The generic form of the Lorentz structures LABXi, that are part of the Feynman 
rule ig2CABXi LABXi for the XaXbX^ three-point vertex in Fig [91 We assume all four- 
momenta, kA,kB,ki, to be outgoing at the vertex. The case of X^X^Xj = VaSV^ is 
trivially related to the case VaVjs S. 



s 



I- (pf ) = ^2«^Ji* ^ ^Zl'Z'J* . (50) 



cos l3Mw 2 



2V2cos/3Mi 



w 



The coupling factors with fS = uD, uL refer to the interaction of an incoming up-type 
quark (u^) or neutralino (z/^) with a Dj- or Lj-sfermion and an outgoing xt ■ In case 
of fS = dU, W, the coupling factors in ( l50l) are related to interactions of an incoming 
down-type quark (d^) or lepton (f) with an Uj- or z/^-sfermion and an outgoing xt^ 
(denoting the charge-conjugate field of xt, see [33]). 

For the specific CABXi factors that emerge at the three-point vertex of the XaXb 
particle pair with the single s-channel exchanged particle Xi in Fig. [HI we refer to the 
Feynman rules in [33]: a specific CABXi is obtained as the factor that multiplies the 
structure ig2LABXi in the respective Feynman rule therein. The generic forms of the 
Lorentz structures Labx, are collected in Tab. |H Finally, (axial-) vector and (pseudo-) 
scalar coupling factors r and q of two SM fermions with a gauge- or Higgs-boson (see 
Fig. nU]) can be directly taken from the corresponding Feynman rules in [33] • 

In order to illustrate how the above rules should be applied, let us consider an exam- 
ple. Suppose we wish to know the coupling factors c^^l^y, a = 1, ... 4, of diagrams with 
s-channel exchange of a Z-boson for the x^jXj^ ~^ W'^G~ XciXea processes. Following 
the recipe above, we draw all tree-level diagrams for the XeiXit, ~^ W~^G~ as well as the 
Xe4,Xe3 ~^ W^G~ process and assign the same fermion flow as given in the corresponding 
generic diagrams. Fig. O Referring to that fixed fermion flow, we determine the follow- 
ing vertex factors in the diagrams ti and s, associated with tree- level XeiXt2 ~^ W~^G~ 
annihilations: 

ym) ^ G , yHti) ^ ^» (_yW* , W* \ 
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= 7" (-<e. + <..75) , CWGZ = (51) 



The coupling factors are those from (jHj) and (H5il . Note that there is no t-channel 
exchange diagram t2 for the above process, as it is forbidden by charge conservation. 
Further, note that the sign in front of the vector-coup hng factor in V'^^lf' and Ve^el 
follows from rule 26) above. In case of diagrams contributing to Xe4Xe3 ~^ W~^G~ we 
find 

^ eAn ''tie4 Fiie4 /5 ! ^ e^ix I Katies, ^ ^iie-i 1^ j ) 

ya{t2) ^ ( W* , W* \ ym) ^ G , G 

Ki? = 7^ {v^Z + 4?e!75) , c^GZ = -f- . (52) 

To obtain the building blocks for the non- vanishing c^") y with a = 1, one has to combine 
the coupling factor expressions in the first row of fIFI]) (the factors related to diagram t\ 
in x^iX^2 ~^ W^G~ annihilations) with the coupling factor expressions in the last row 
of (!52|) (expressions originating from diagram s in Xe4Xe3 ~^ W^G~^^ as a = 1 refers to 
the t\s product of tree-level diagrams. Similarly, for a = 2 and 4, one has to build the 
combinations of expressions referring to st\ and 5^2- Therefore, the building-blocks for 
the non-vanishing c>^\^y related to single s-channel V = Z exchange read: 

« = 1 : {{^n, , {-C:^ <.} ^ {«f ^ ^ {-4^}} ^ (53) 
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« = 4: {{-<e„<ej, {s%,.V%,}^ {-^4,^4}' {-$}} • (55) 

In selecting one element from each of the above given subsets and multiplying the selected 
elements for fixed a with each other, the Cn\v expressions in x^^x[^ W'^G~ — )■ Xe4Xe3 
reactions are obtained. Proceeding in that way, we obtain eight different coupling factor 
combinations for fixed a, that are labeled with index n. Following our convention for 
this label, n ranges over n = rrr, rrq, rqr, qrr, rqq, qrq, qqr, qqq. The Cgg|, ^_^y expression, 
for example, reads 

(2) _ _Sw_ Z W G* /c-f^N 

'-qqr,hV " ^2 "eiea "iieg ■^1164 " 

c — 



w 



A. 3 Kinematic factors 

The kinematic factor expressions that refer to a specific XeiXe2 ~^ ^a^b Xe^Xes 
scattering reaction depend on the external particles' mass scales m, m, and M = m + m. 
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We remind the reader of our convention (see Sec. 12.41) 



rriei + rrie^ _ irie^ + meg 
m = , m = , 

2 ' 2 ^ ^ 

Further recall that we expand the scattering amplitudes in dm, 6fn and count these 
quantities as (!?(w^). Hence, for the leading-order S'-wave results presented below, the 
mass differences 6m = 5m = 0, such that there are only two mass scales, m and m, 
left, which characterize the external chargino or neutralino states. The masses of the 
particles Xa and Xb will be denoted with and m^. Let us introduce the general 
notation for the rescaling of any mass in units of the mass scale M, 

= ^ . (58) 

Define the dimensionless quantities 

^AB = m\-m\ , 



(3 = ^l-2{m\ + ml) + A\^ , (59) 

where /3 is the leading-order term in the expansion of the velocity of the XaXb particle 
pair in the non-relativistic momenta and mass differences. The expansion of single 
s-channel (gauge or Higgs boson Xi) exchange propagators in 5m, 5m and the non- 
relativistic 3-momenta leads to the following denominator-structure at leading order: 

Pi = 1-^? • (60) 

Similarly, the leading-order expansion of t- and w-channel gaugino and sfermion prop- 
agators in 5m, 5m and the non-relativistic 3-momenta gives rise to the denominator- 
structures 

Pi,AB = m m + m^ — m m\ — m m\ ; 

Pi.BA = Pi,AB \a^b ■ (61) 

Using the above definitions, the kinematic factors for the leading order S'-wave Wilson 
coefficients related to the selfenergy-topology in Fig. |5] are conveniently written as 



B^,t^,-C'^'Sj) = ^ B^,t^fC'^'Sj) , (62) 



where the labels ii and 22 refer to the particle species that are exchanged in the left and 
right s-channel propagator. As generically either gauge-boson {V) or Higgs (S) s-channel 
exchange occurs in the processes under consideration, the combination iii2 is given by 
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^1^2 = W, yS, SV, SS. Kinematic factors arising from the triangle-topologies shown in 
Fig. [6] have the following generic form 



Ci^W, (2.+!^^) ^ Ci^aW, (2.+!^^) ^^3^4 (g3) 

The index ii in the above expressions is related to the t- or u-channel exchanged par- 
ticle species, whereas the subscript-index X indicates the type of the single s-channel 
exchanged particle-species, X = V, S. Finally, kinematic factors associated with the 
box-topologies of Fig. [7] and Fig. |8] generically read 

t^{1)XaXb f2s+l c \ _ XaXb ^2s+l Q \ 

j^{2)XaXb (2s+l n \ _ n(2)XAXs/2s+l Q \ 

-< iiAB ^i2AB 

j^{^)XaXb (2s+l n \ _ 7=)(3)XAXfl /2s+l Q \ 



Indices ii and 22 in refer to the exchanged particle species in the left and right t- 
and u-channels of the 1-loop box-amplitudes, respectively. 

Throughout this appendix, the labels A and B are related to the particles Xa and 
Xb- Recall that these are the actual final-state particles in a XiXj ~^ ^a^b (tree- 
level) annihilation reaction. The overall prefactors in f l62] - [64|) arise from the phase- 
space integration (/3) and from the leading-order expansion of s- or t- and u-channel 
propagators in the non-relativistic limit. 

Finally recall, that each individual index n in fl52HM|) is given by a character string, 
whose elements indicate the type (r or q) of the corresponding generic coupling structures 
at the vertices of the respective underlying 1-loop amplitude in Figs. [SHHl In the results 
that we quote next we only write explicitly the kinematic factors for those n which are 
non- vanishing. 

A. 3.1 Kinematic factors for XaXb = VV 

The kinematic factors B^\i^ in case of ^Sq partial wave reactions are given by 

B^.^vvCS,) = - ^ + 3 Ai^ , (65) 
B^J,vs{'S,) = B^.'^svCSo) = 3 fhwAAB , (66) 
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B^^CSo) =^ml. (67) 
In case of ^Si partial-wave reactions we have 

= - l^\P"-\^"^B- (68) 

Only the kinematic factors B^\i^ given explicitly in (pHH^Sj) with n = rr, qq are non- 
vanishing. In case of X^Xb = VV, the kinematic factors for the triangle- and box- 
diagram topologies a = 3,4 are related to the corresponding expressions for diagram- 
topologies a = 1,2 (see Figs. [6H7]). The relations read 

^{3)VV /2s+l c \ _ r^WVV f2s+l Q \ I 
^{4)VV f2s+l Q \ _ /=((2) VV/2S+1 Q \ I 
-=y(3) yy /2s+l c ^ _ r^WW /2s+l n \ I 

/=y(4) yy /2s+i Q \ _ /=y(2) yy /2s+i c n i 

^n,hi2 I ^J) - ^n,ni2 \ ^J) \A^B , 

7=)(4) yy /2s+i Q \ _ n(2) yy /2s+i r. \ i ff^Q\ 

The minus sign in the relation for the triangle coefficients Cl^l\^ in (169|) arises from 
interchanging the two gauge-bosons Xa_ and at the internal three-gauge boson vertex. 
The expressions C"^"- diagram-topologies a = 1,2, that refer to leading-order ^5*0 
partial waves read 

?2 



(i)yy /ic, X _ /3 3 ,^ 3 



C'Z:M'So) = y - I (m - m)AAB - | + 3 m,, A^^ , (70) 
^(2) yy /I X _ ^(1) yy .1 r, N z^-, x 



In case of C„ expressions related to '^Si partial waves and diagram-topologies a = 1,2 
we find 

Ci^lZCS,) = - ^ + (m - ^)^ + ^ + 3 . (72) 



We deduce the following expressions for C^l coefficients and diagram topologies 
1,2: 

Cr^Jg^i^sC So) = CgfZ^sCSo) = -mw{m -m + Aab) + 4 mwfhi, . (73) 

There are additional non- vanishing C^^]^^ expressions for X = V, S in both the case of 
and ^S'l partial wave reactions, which are related to the expressions ([70]) and fl72]) in 
the following way: 

^(i)yy /2s+1q \ _ ^(2)yy /2s+i c n _ M'i-)vv /2s+1q \i 



a 
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Ml)VV /2S+1Q \ _ M'i)VV /2S+1Q \ _ r'WVV /2s+1q \\ (rj.\ 
^qrq,hX\ J ) " '-'qrq,iiX{ J ) " '-'rqq,iiX{ ^ J ) \mi^-^-m,^ ■ [I'i) 

Turning to terms related to box- diagrams, the non-vanishing expressions 

^iir for 

a = 1,2 are given by 



rhi^ {m- m- Aab) - fhi^{fh - m + Aab) , (76) 
- {fhi^ + fhi^) {m -m + Aab) , (77) 



and 



= - ^ - ^ (m - M)2 + ^ A^^ + 2 m,,m,, , (78) 

^il'^LC'^i) = l^'-l (m-m+ Aab)' - 2 m.^m.^ , (79) 
Dizzy's,) = - 1 - ^ (m - M)2 + i Ai^ - 2 m,m, 

+ mi^{m-m- Aab) +'mi^{'m-m + Aab) , (80) 
^i'i^'Le^i) = -^/3' + ^(m-m + A^B)' + 2m,m, 

- mi^{m -m + Aab) - mi^{rh -m + Aab) ■ (81) 

The remaining non-vanishing kinematic factors D^\i^ for both spin-0 and spin-1 XX 
states are related to the expressions given above by 



qqqq, tli2 


[^^^^^.) 


_ fjia)VV 

rrrr, i\i2 


C'-^'Sj) 


5 


rrqq, 


('^+^^j) 


_ f^ia)VV 
qqrr, 1112 


C'^-'Sj) 


_ f)io')VV /'2s+lr' 
rrrr,iii2y J 




(^^+^^i) 


_ p,ia)VV 

rqqr, 11*2 


C'^-'Sj) 


1 "lii,2~^ ~™'n,2 ' 


rqrq, 


[^'^'Sj) 


_ f)ia)VV 

rqqr, i\i2 


C'^'Sj) 


1 rhi^-^ -™i2 ' 




C'^'Sj) 


_ fjia)VV 
rqqr, i\i2 


C'^-'Sj) 


1 rhi-^-^ -rhi-^ ■ 



(82) 

The notation in the second line of fl82|) means that the product fhi-^rhi^ is replaced, but 
all other appearances of either or fhi^ are untouched. 
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A. 3. 2 Kinematic factors for XaXb = VS 

We find the following expressions for -B^f^j^ terms in ^Sq partial- wave reactions with 
t^t^ = VV, VS, SV, SS: 

B^'vvCSo) = -ml^ (83) 

Bl'vsi'S,) = Bl%yCSo) = ^ (-3 + A^^) , (84) 

^S.(^^o) =^-1 + IAas-^. (85) 
In case of ^Si partial-wave processes the corresponding -B^f^jj coefficients read 

Bl^yCSi) = ml, . (86) 

Kinematic factors C^l \f , that are related to the four generic triangle-topologies a with 
gauge-boson [V) exchange in the single s-channel (see Fig. |6]) are given by 

ClZU'So) = Cr^C^o) = -%^(m -m+AAB) - mwm,, , (87) 

ClZU'So) = CSC^o) = "^{fh-m- Aab) + mwm. . (88) 

as well as 

/=<{3) /3 r. N _ /=y{4) VS /3 Q N _ r^(3) VS /l g ^ , ^ /q^n 

^rrr,iiV\ ~ ^rrr,iiV\ '-'^J ~ ^rqq,iiV\ ^0 J | m^^ -s.-mi^ • l^^U; 

In case of C^\^g expressions we find 

Crqq%{^So) = - ^ - ^ (fh - W) + {fh - W - 3) 



■AB 



+ %-%(3-A^^), (91) 

^qqr,iiS\ >^0) - ^rqq,iiS\ '^Oj , (.y^j 

+ % + %(3-A^^), (93) 

C^frZCSo) = ClZ'sCSo) . (94) 
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There are additional non- vanishing kinematic factors for C^l^x "with X = V or S, 
related to the corresponding expressions in fl57HM|) in the following way: 





'2s^ 




^i2)VS 
- ^rqq,hX^ 


'2sH 




_ _ Mi)vs . 

rrr, iiX'' 


■2s + lq \\ 

J J )\mi-^^-mi^ 


M3)VS 
^qqr,hX^ 


'2sH 


-'Sj) 


_ Mi)vs 

~ '-^rqq,iiX^ 


'2sH 


-'Sj) 


_ M3)VS /2s- 
rrr, iiX\ 




^qrq,hX^ 


'2sH 


-'Sj) 


_ M-i)vs 

~ ^qrq,iiX^ 


'2sH 


-'Sj) 


_ Mi)vs 

~ '^rqq,iiX^ 


'2s+l c \\ 


^qrq,hX^ 


'2sH 




^(4)ys 

~ ^qrq,iiX^ 


'2sH 




^{3)VS ,2s- 
~ ^rqq,hX\ 


Sj)\fhi^-^-fh.i^ ■ 



(95) 

The non-vanishing kinematic factors for Xj^Xb = VS and the four box-topologies a are 
given by 

+ ^rhi^{m-m- Aab) + ]^ mi^{m - m ^ A^s) , (96) 



-]^{mi^+mi^){fh-m + /^ab) ■ (97) 

-^rgg^fl 43(^*^0) = -D^gg^fii2(^'S'o)UoB , (98) 
-^rgg^fi 12(^*^0) = Df^^^^^^^^i} Sq)\a^B ■ (99) 



In case of ^5*1 partial waves we have 

— - - (m - mf + - 
12' 4 ^ ^4 



- ^ rhi^ {m-m- Aab) + ^ fhi^ifh - m + Aab) , (100) 

+ ^ {mi^+fhi^){m-m + Aab) , (101) 

rrrr, 1112^ rrrr, 1112^ l/lm-<->m ' V ^ / 

r)(4)l^S .Son ^ r)(2)V5 /3 o N 1 ^ /inoN 

rrrr, jii2V "-^ly -'-^ rrrr,i\i2\ ' V ^"-^Z 

Relations for the remaining kinematic factors for both ^5*0 and '^S'l read in case of 
diagram-topologies a = 1,2: 

f,(a)VS /2s+lq \ _ (_T\a f^{a)VS /2s+l r- \| 

^qqqq,iii2\ J I ~ \ ^) -^rrrr,iii2\ '-'Jjl '^ii,2^ ~^h.2 ' 
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f){a)VS /2s+1q \ _ (_T\a+l f){a)VS /2s+l r. \\ 

f){oi)VS f2s+ln \ _ _ f){o')yS (2s+lq \\ 

^qqrr,iii2\ "^JJ ~ ^rrrr,iii2\ '~'J) \ rrii-^-^ -nii^ ; 

f)(a)V5 /2s+lc N _ (_T\a f,ia)VS /2s+l C \\ 

^qrrq,iii2\ J J \ ^) ^rqqr,iii2\ '-'Jjl ^ii^2^ ~™n,2 ' 

f){a)VS /2S+10 \ _ fS^")^"^ (2s+l o \| 

^rqrq,iii2\ '-'J) \ ^) ^rqqr,iii2\ "^Jjl mi^-^ -''"i2 ' 

f\{a)VS /2s+1q \ _ _ n{")^'S' |'2s+lc< \| 

^qrqr,hi2\ '~'Jl ~ ^rqqr,hi2\ '~'J)\ m^.^^ -rhi^ ■ \^^^ ) 

The corresponding relations for diagram-topologies a = 3,4 are given by 

f^{a)VS /2s+lq \ _ ( f){a)VS /2s+l r. \ \ 

^qqqq,ili2\ "^JJ ~ \ ^) ^rrrrAii2\ '^J)\ ™»i,2~^ ~™H,2 ' 

p.{a)VS /2s+lq \ _ f){a)VS /2s+l C \\ 

^rrqq,iii2\ '-'J) \ ^) ^rrrrAii2\ '-'J)\ nii^^ -"^^2 ' 

f\{a)VS /2s+lq \ _ jj{a)VS /2s+l q \i 

^qqrr,iii2\ "^JJ ^rrrr,ni2\ '~'J)\ rhi^^ -rhi^ ; 

f)ia)VS /2s+lr. \ _ (T \a f)ia)VS /2s+l c \\ 

^qrrq,iii2\ J > ~ \ ^) ^rqqr,iii2\ '^J)\ "^ii,2~^ ~™n,2 ' 

){a)VS (2s+lq \ _ ( iNa f\{a)VS /2s+l , 



^rqrq,hi2\ '-'J) ~ \ ^rqqr,hi2\ '^J)\ rn 

f){a)VS f2s+lq \ _ p,ia)VS /2s+l r. \\ H flfi^l 

^qrqr,iii2\ '-'J) ~ ^rqqr,iii2\ '^J)\ mij^ -mii " [lUO) 

A. 3. 3 Kinematic factors for XaXb = SS 
The non-vanishing B^^^_^-^ terms with ziZ2 = VV, VS, SV, VV read 



fySS 1 
^qq,VV^ 


['So) 


- 


tdSS I 
^qq,VS^ 


:'so) 


= Bqq^SvC^o) = — fnw'^AB 




'}So) 


= fn^w ' 



and in case of ^Si reactions 



(107) 
(108) 
(109) 



B'JM'Si) = y • (110) 

As in the case of X^Xb = VV, the kinematic factors for X^Xb = SS and diagram 
topologies a = 3, 4 are related to the corresponding expressions that arise from diagram- 
topologies a = 1,2. This applies to both triangle- and box-topologies (see Fig. [6] and 
Fig.m): 

/=((3) 55/2S+1 c N _ ^{l) SS f2s+l q 



-(4)S5/2s+lQ N _ r'{2)SS(2s+lq \ i 
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) U<->B 


7 






) \a^b 




fS{3) 5S/2S+1 o N 


A(l)5S/2s+l o 


\ 1 

j \a^b 


5 


7=)(4)5S/2s+1q \ 


A(2)SS/2s+lr- 
--'^n, ilia I '^'^ 


) 




Cl^lv^ for dia^ 


rram-topologies a 


= 1,2 


we find 


_ /V(2) SS fl Q \ 


= m — m 

9 ^ 


+ A^B 


) + rhi^ 



^S,1fv('^o) =C:^;'rT.vCSo) = ^{m-rn + A^B)+mAAB, (112) 

The C^^l^g^ expressions with a = 1, 2 are given by 

CiZU'So) = C'^^!S('^o) = (m - M + A^^) - m^m, . (114) 

All other non- vanishing expressions for Ci^'lf^ with X = V,S and a = 1, 2 can be related 
to the terms in flll2Hll41) in the following way: 

Mi)SS f2s + lq \ _ ^{2)SS f2s + lq \ _ _ Mi)SS f2s + lq \\ 

'^qqr,iiX\ '-'J J ~ '^rqq,iiX\ '-'J J ~ '^rrr,iiX\ ^ Jj I m^^ -5>-m,-^ , 

M^)SS ,2s+lq \ _ ^{2)SS /2s+1q \ _ ^WSS ,2s+l q \\ 

'~'qrq,iiX\ J ) — '~'qrq,iiX\ J ) — '~'rqq,iiX\ J J ) \mi^^-m,-^ ■ l-LJ-OJ 

The expressions D'^n^it diagram-topologies a = 1, 2 and ^5*0 partial waves read 

+ ^ (^-^- Aab) + ^ (m-^+ A^b) , (116) 



+ ]^{mi^+mi2){m-m + Aab) ■ (117) 



In case of a ^5*1 partial wave configuration we find 



2 



'^{a)SS (Zq \ _ ( 1^Q ^ 

12 



DTJ:X.,{'S,) ={-\Y'-. (118) 



The remaining non- vanishing kinematic factors related to both ^5*0 and partial-wave 
reactions read 

f^{a)SS f2s+lq \ ^ f){a)SS /2s+lr< N 
qqqq,ili2\ '-'J) rrrr,iii2\ ^J) : 
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jjia)SS 
rrqq, iii2 


f2s^ 


~'Sj) 


_ jj{a)SS /2s+lr. \ _ 
qqrT,i\i2\ '-'J J 


- -D 


ia)SS /2s+lr' \ 
rrrr,iii2\ J J 1 




f\ia)SS 
qrrq, 1112 


f2s^ 


-'Sj) 


_ jjia)SS /2s+lr. \ 1 
rqqT,i\i2\ J) ' 1 








T\(a)SS 
rqrq, 1112 


f2s^ 


~'Sj) 


— _ jji'^)SS /2S+10 N 
rqqr, 1112^ J J 


1 rhi^- 






jj[a)SS 
qrqr, 1112 


/2s^ 




— — f)('^)^^ (2s + lq \ 

rqqr,iii2\ ^J) 




* ~rhi^ ■ 


(119) 



A. 3. 4 Kinematic factors for XaXb = ff 

The non-vanishing B^/^^^^ terms with iii2 = VV, VS, SV, SS are given by 

BiLvvCSo) =1-/3^ + 4 m^ms - , (120) 
Bgq,vsCSo) = Bg^,svCSo) = 2 (m^ + - (m^ - m5)A^s) , (121) 
Bgq,ssCSo) = 1 + + 4 inAfhB - A\s , (122) 

and in case of ^5*1 partial-wave reactions 

HUvvCSi) = 1 + y + 4 ruAmB - A^^ . (123) 

There are additional non- vanishing terms -B^'^i^j^ I'slcited to the expressions in fll20H123"]l . 
In case of «i«2 = VV, SS, the corresponding relations read 

P>ff (2s+lq \ _ hff f2s+lci \\ _ _ _ ri9zll 
^rqqr,ni2\ J ) ~ ^rrrr,iii2\ ^ J ) ImAniB^-mAmB ) K^^^^J 

B(rrq,iii2( ^j) ~ Bgqgq^i^i^i S j)\ffi^ffig^-fhArhB , (125) 

and our notation implies, that the product mArriB has to be replaced, but all other 
occurrences of or rriB are untouched. Similarly, in case of iii2 = VS,SV, the 

~ f f 

additional non-vanishing BJ^^^^ terms are given by 

Brqqr,iii2( ~^ ^j) ~ ~ -^r/rr, 1112 ( ^j)\mA^-mA 5 (126) 
Bqrrq,ni2{ ^ ^ j) ~ ~ B^qq ^^^^{ S j)\fnA-^-mA ■ (127) 

Here our notation means, that niA is replaced, but the friA in ^ab is left unchanged. The 
relations among kinematic factors for diagram topologies a = 3, 4 and diagram-topologies 
a = 1, 2 in both the cases of box- and triangle-topologies are given by (X = V, S) 



^(3) // , 




^(1)//, 






(128) 


^(4) // , 


'^'^^Sj) 


^(2)//, 


"'^'Sj) 


\a^b , 
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n, i\i2 ' 


:''^'sj) 


n, i\i2 


C'^'Sj) 




(130) 
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D 



(4) fff2s+l 



n, 1-112 ^ 



n, ii%2 ^ 



(131) 



compare to the generic diagrams in Fig.[6]and Fig. El The structures Cl^lfy for topologies 
a = 1,2 are given by 



^qqqq,ilV\ '^O) 



y-^ (l-2m^)(l-2mB) 



{niA - rriB) 



A 



AB 

4 ' 



and in case of partial wave reactions the respective expressions read 



^-^ (l + 2m^)(l + 2ms) 



+ {friA - friB) 



A 



AB 



^AB 



(132) 



(133) 



The relations of the additional non-vanishing c'^\^y expressions to the respective terms 
in f lT32HT33|) read 



/=.(!)// 1 

qqrr, ii V ' 


'2s^ 


-'Sj) 


_/=^(2)// , 
rrqq, iiV' 


'2s- 
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' '^qqqq,hV\ J > 
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The terms C^l'^J for a = 1, 2 read 



qqqq,hS\ 



So) 



^ + \ (l-2m^)(l-2mB) 



+ {friA - ms) 



A 



AB 



^AB 



(135) 



and all remaining non-vanishing C^"j^// terms are obtained from f ll35p in the following 
way: 



^rqrq,nS\ '~'J> 
''qrrq, iiS\ 



/=<(2)// /2s+1q \ _ M»)ff f2s+lq \\ 



''qrqr, iiS\ 

=<(")// / 
^qqqq,ilS\ 
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rrqq,iiS\ '-'■J J ~ ^qqrr,hS\ J ) ~ ^qqqqfhsi ~^ ^ j)\mA,B^-mA,B ■ (136) 



In case of box-diagram topologies a = 1,2, we find the following D^^]^/^ structures for 
the ^5*0 partial waves: 



rrrr, i\%2 * 



(137) 
(138) 



For ^Si partial-wave configurations we have 



rrrr i\i2 



Si] 



rrqq,iii2^ 1 



rrrr, iii2 ^ 
?2 



^o) 



(3 



-1)"+M ^ + - m^ms . 
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(139) 
(140) 



Relations for the remaining non-vanishing kinematic factors related to both ^5*0 and ^Si 
partial-wave processes read in case of diagram topology a = 1 read 



^qqqq, ii?2 




= £)«^^ . 1 

rrrr, i\i2 




rhA,B^ -rhA,B ' 


D^'^". . , 

qqrr, i\i2 




rrqq, 1112 






rqqr, i\i2 




= b^')ff. . 

qrrq, 1112 




— f)^^^^^ (2s + ln 
rrqq,iii2\ J 


b^'^ff. . 1 
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1 m^— > — "lA ' 


qrqr, i\i2 




= D^'^" . 1 

rrrr, i\i2 




rhs-^ -rriB ■ 



-mATriB 1 
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In case of diagram topology a = 2, the corresponding relations are given by 



qqqq,i\i2 ^ 



qqrr, Jii2 * 
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D 



(2)ff (2s+l 
rqqr, 1112 \ 

qrrq, iii2 V 



Sj) 



s. 



rqrq, i\i2 ^ 



rrrr, i\i2 ' 



J)\ mA,B->- -rh,A,B ' 



D 



(2)// /2s+l 
rrqq, 1112 \ 



s, 



n(2)// (2s+lq \\ 

I ^J)\ ifiA^ —rriA 1 



rrrr, i\i2 ^ 



D 



mff /2s+lr- 

rrrr,iii2\ '-'J)\ rhs^ -niB ) 



qrqr, i\i2 ' 



D 



(2)ff (2s+l 
rrrr, 1112 v 



Sj) 



mAiriB—^ —niArriB 1 



(142) 



A. 3. 5 Kinematic factors for XaXb = rj-q 

In case of XaXb = Tjfj one cannot directly construct the coupling factors 6„ i^i^ using 
the recipe given in Sec. IA.2t which is based on considering the XeiXe2 ~^ XaXb and 
Xe4,Xe3 XaXb tree-level annihilation amplitudes. In order to obtain the coupling 
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factor expressions bniii2, that correspond to the kinematic factors presented below, one 
should proceed as follows: First extract the (axial-) vector and (pseudo-) scalar coupling 
factors associated with the interaction of the XeiXe2 or Xe^Xe-j pair and the s-channel 
exchanged particle species. This is done following the steps 1. and 2. in the recipe 
given in Sec. IA.2I Next, complex-conjugate the couplings related to the Xe4Xe3 particle 
pair. In order to determine the couplings to the ghosts, consider the 1-loop amplitude 
XeiXe2 VV ^ XeiXe-i, similar to the selfenergy- amplitude in Fig. |5l Assign a ghost 
flow to the lower line of the 1-loop amplitude (labelled with in Fig. [5]), that flows 
from left to right. Consequently there is a ghost flow from right to left on the upper 
line, which is labelled with Xb- Assume that the coupling factors at each of the two 
ghost vertices are generically of the form ig2CABXiLABXi, where the Lorentz structures 
LABXi are defined in Tab. HI Determine the expressions that replace the generic CAsXi 
factors for the specific process under consideration. Now build all possible combinations 
of two-coupling factor products from the set of the neutralino/chargino couplings to the 
s-channel exchanged particles (including factors of —1 in front of vector couplings) in 
the XeiXe2 Wi ^ XeiXe^ reactlou, and multiply them by the CABXi^ and CAsXi^ factors. 
The convention for the naming of the resulting coupling factor expressions hn,i^i2 with 
subscripts n = rr, qq is the same as in the cases XaXb = VV, VS, SS, see Sec. IA.2I 
The coupling factors j^jj derived in this way correspond to the kinematic factors given 
below. Note that the mass parameter ruA in the expressions below refers to the mass of 
the ghost flowing in the lower line, and ms to the mass of the ghost in the upper line. 
The non- vanishing -B^^j^ terms with iii2 = VV, VS, SV, SS read 

BZvvCSo) =\{l-AU, (143) 

B^vsCSo) = -%^(1 + A^b), (144) 

BZsvCSo) =%(1-A^b), (145) 



^qg,SV\ 



BZssCSo) = -ml,. (146) 



Similarly, 



Bl^M'Si) = - (147) 
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